ERROR TERM IMPROVEMENTS FOR VAN DER CORPUT 

TRANSFORMS. 



J. VANDEHEY 



Abstract. We improve the error term in the van der Corput transform for expo- 
nential sums 

. — , * 

g(n)exp(2mf(n)). 



a<n<b 



For many functions g and /, we can extract the next term in the asymptotic, 
showing that previous results, such as those of Karatsuba and Korolev, are sharp. 
Of particular note, the methods of this paper avoid the use of the truncated Poisson 
formula, and thus can be applied to much longer intervals [a, b] with far better 
results. We provide a detailed analysis of the error term in the case g{x) = 1 and 
f(x) = (x/3) 3 / 2 . 

1. Introduction 

We are interested here in estimating the error term A associated with the van der 
Corput transform, 

m V* / \ (fl u g(x r )e(f(x r )-rx r + \) 

(1) 2^ g(n)e(f{n)) = ^ 7WT7 ' 

a<n<b f'(a)<r<f'(b) V * ^ Xr ' 

where / and g are several times continuously differentiable functions with f"(x) > 
for x G [a, b] and where x r is defined by f'(x r ) = r. A starred sum indicates that if 
a limit of summation is an integer, the corresponding summand is multiplied by 1/2. 
The function e(x) denotes e 2mx . 

The van der Corput transform is best-known in number theory for being the cru- 
cial element of Process B in the theory of exponent pairs and is sometimes simply 
referred to as "Process B." As a method of estimating exponential sums, the van 
der Corput transform is often presented alongside other methods of Weyl, van der 
Corput, and Vinogradov. Direct application of the van der Corput transform can 
take a complicated sum to one more amenable to estimation techniques or it can re- 
duce the number of terms and therefore make computational estimations easier. The 
van der Corput transform itself is involutive — applying it to the right-hand side of 
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([I]) will simply return the sum on the left-hand side of ((!]) — so one gains nothing by 
applying it twice in a row; but one could alternate applications of the van der Corput 
transform with other techniques (such as the Process A of the theory of exponent 
pairs) to achieve better results. This alternation method is still a fruitful ground 
for modern research. Recently, Cellarosi [2] attained interesting new results in the 
classical case of g(x) = 1 and f(x) = ax 2 , where the alternating technique employed 
is simply reducing a modulo 1; Nakai [2H1 EH I2H] has investigated the possibility 
of using an analogous method when f(x) is cubic or quartic. The van der Corput 
transform has also seen recent use in physical applications, including quantum optics 
and wave processes (see [15] and the papers cited there for more details). 
Van der Corput [33 J originally showed that, given 

|f(x)|xA 2 , |/( 3 )(x)|«A 3 , g(x) = l, for are [a, 6], 

the error term can be bounded like 

A = 0{\~ 1/2 ) + 0(log(/'(6) - f{a) + 2)) + 0((b - a)Aj /5 Af ). 

(Here we use the Landau and Vinodradov asymptotic notations which will be de- 
fined explicitly in Section [2]) Phillips |29j improved this error under the additional 
assumptions 

|/ (4) (x)| < A 4 for x G [a, b] and A3 = A 2 A 4 ; 

in this case, we can replace 0((b — t^A^A^ 5 ) with 0((b — a)X^ 3 ). 

The form of the error term found in most modern texts on analytic number theory 
[131 EHl E3fl has its roots in the work of Kolesnik |19j and Heath-Brown [ID] , 
although the results of the latter authors required analyticity of the function /, an 
assumption which has since been circumvented. This moderately-difficult-to-prove 
form of the error term suffices for many basic application of the van der Corput 
transform. We present this modern bound on the error in the notation used by 
Huxley. 

Theorem 1.1. (Lemma 5.5.3 in [12] ) 

Suppose that f(x) is real and four times continuously differentiable on [a,b]. Sup- 
pose that there are positive parameters M and T, with M > b — a, such that, for 
x G [a, b], we have 

f"(x) x T/M 2 , f 3 \x) < T/M 3 , and f {A \x) < T/M 4 . 

Let g(x) be a real function of bounded variation V on the closed interval [a,b]. Then 

g(x r )e(f(x r ) - rx r + |) 

a<n<b f'(a)<r<f'(b) 



E «K»M/(»» = ^ r 



1 Curiously, [35] skips this form of the error term entirely. 
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+ O [ (V + \g(a)\) [ *L + l g (/ '(6) _ f {a) + 2 ) 



where x r is the unique solution in [a,b] to f'(x r ) = r. The implicit constant in the 
big-0 term depends on the implicit constants in the relations between T, M , and the 
derivatives of f(x) . 

The error term Mj \pT here corresponds to the A 2 ^ 2 term in the estimates of van 
der Corput and Phillips. 

Unfortunately, for many interesting cases, the above error is insufficient. As Huxley 
[121 p. 475] notes, when applying the van der Corput transform to a multi-dimensional 
exponential sum, "...the error terms and the truncation error in the Poisson summa- 
tion formula may add up to more than the estimate for the reflected sum." Thus, 
finer error terms, useful for a broad spectrum of problems including computation and 
physical applications, have been given by various people, including Kolesnik [IB], Liu 
[21] . Redouaby and Sargos [31], and Karatsuba and Korolev [15] . Liu extends an 
(unfortunately obscure) earlier work of Min [22] . removing the latter's condition that 
f(x) be an algebraic function. Redouaby and Sargos show that the conditions on 
f^(x) and g"(x) can be removed without greatly increasing the bound on the error 
term. The work of Karatsuba and Korolev is unique among all papers cited here, in 
that it is the only to give the implicit constants in the big-0 terms explicitly, making 
it superior for most computational work. 

Outside of Redouaby and Sargos' result, we cannot briefly state any of these other 
forms of the error term in full detail; we will, however provide the following inexplicit 
form of Karatsuba and Korolev's result as an example of the comparative strength 
of these errors compared with Theorem 

Theorem 1.2. Suppose that f(x) and g(x) are real-valued functions with f four 
times continuously differentiable and g two times continuously differentiate on the 
interval [a,b\. Suppose there are positive constants M, T, and U, with M x b — a, 
such that, for x £ [a, b], 

(2) f"(x) x T/M 2 , |/( 3 )(a;)| « T/M 3 , \f^(x)\ « T/M 4 

(3) \g(x)\ « U, \g'(x)\ « U/M, \g"(x)\ « U/M 2 . 

Then, 

E*»(»W(»))= E* g(l - 0e(/(l ;':7 + 1/8) +o(t/(r( a ) + T W) ) 

a<n<b f'{a)<r<f'{b) V ■> \ Xr ' 

+ o(u (log(/'(6) - f\a) + 2) + f + ^ + 1 
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where x r is the unique solution to f'(x r ) = r in the interval [a,b] and 

o, \\f'(ri\\ = o 

II/'(aOII^o' 



(4) T(ji) = . fM 1 

1 mm ' 



t wmw 

The size of the implicit constant in the big-0 term depends on the implicit constants 
in the relations of M, T, U , and the derivatives of f(x) and g(x). 

If, in addition, M < T < M 2 and M > b—a, we may remove the terms 0(T/M 2 + 

1)0 

The M/y/T term again makes an appearance in this theorem. It cannot, in fact, 
be completely removed. If b is not an integer, but f'(b) is an integer, then changing 
b to b — e for some very small e > will not change the value of the sum on the 
left-hand side of ([T|, but removes a term of size M/yT from the right-hand side. 

Van der Corput's results in [31] also deserve a mention here, as they are of an 
entirely different flavor from those listed above and because they are generally not 
very well-known^ Instead of having a coefficient of g(x r ) / \/ f"(x r ) in the transformed 
sum, he has a more general set of coefficients. In addition, his error term — which 
would require too many definitions to state succinctly here — bears no resemblance to 
any of the other error estimates cited or formulated in this paper. 

While van der Corput's results in [31] are quite complicated to use, they are aesthet- 
ically pleasing. As we remarked above, the van der Corput transform is involutive, 
but in all the other results given above as well as the main result of this paper, one 
can obtain very different error terms when one applies the transform to the right- 
hand side of ([T| instead of the left, assuming that the conditions necessary to apply 
the results would even still hold. In [31], van der Corput shows that the conditions 
are still satisfied and the error term unchanged regardless of which side one applies 
his transform to; his theorem is involutive. 

The van der Corput transform and its error has been studied in much more general 
settings than we go into here: of particular interest, Jutila [H] considered sums of 
the form ^b(n)g(n)e(f(n)) for certain multiplicative functions b(n) (see also p2J 
Ch. 20]), and Kratzel [20] considered the van der Corput transform of a convergent 
infinite series. 

One may ask what the best possible error term could be. Given the frequent 
restriction in theorems on the van der Corput transform that f"(x) x A2 (or, equiva- 
lently f"(x) x T/M 2 ), it is not surprising that the case where g is constant and / is 



2 This supplementary result comes from Karatsuba and Voronin [17] . 

3 This may be largely van der Corput's fault. Robert Schmidt, in reviewing van der Corput's 
paper for Zentralblatt, remarked that because van der Corput provided no context for his results, 
neither in terms of past results or future goals, that the paper was unlikely to spark much interest, 
and indeed it has so far only been cited once elsewhere, in [3]. Due to the desire to avoid repeating 
van der Corput's mistake, the complicated main theorem of this paper will not be presented in the 
introduction. 
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quadratic (that is, f"(x) = A2) is the most-commonly studied special case of the van 
der Corput transform and the one with the best error terms [21 El El |25j EH] . Fedotov 
and Klopp [8^] have given the error term in this case as an explicit integral. But 
perhaps most amazing are the results of Coutsias and Kazarinoff [3]: they showed 
that for positive integers n, we have 



N 

E* 

fe=0 



00 



2" 



e(sgn(uj)/8) 



E* 

fc=0 



1 

CO 



2 



< C 



N 



for < \oj\ < 1, N — \n/oJ\ is the nearest integer to n/u, and 1 < C < 3.14 is a 
particular constant. Not only is the error bounded but it shrinks to zero as n/u nears 
an integer. 

The Coutsias-Kazarinoff result suggests that the van der Corput transform should 
be very accurate; in particular, the van der Corput transform for nice enough func- 
tions / and g shouldn't have compounding error terms (such as the log(/'(6) — f'(a) + 
2)) seen in all the other results mentioned above. 



The main theorems of this paper (and Theorem 4.3, especially) confirm this hy- 
pothesis, allowing the van der Corput transform to be applied on very long intervals 
with a much higher degree of accuracy than in previous results. 

As a quick example of this, consider the following well-known transform which 
appears in Iwaniec and Kowalski's book [131 p. 211] (a version of this transform also 
appears in [2T]). Given X > 0, N > 0, and a > 1, v > 1, consider 



(5) E 

N<n<uN 



X 

a 



where l/a + 1/(3 = 1, 
and Kowalskjj give 



n 
N 



■ E* 

M<m<iiM 



1 



J Km) 1 a 



and MN = X. Using a form of Theorem 1.1 , Iwaniec 



A = 0{N- 1 ' 2 log(M + 2) + M" 1/2 log(iV + 2)), 

with an implicit constant dependent on a and v. 

Using the results of this paper, we may improve this to the following. 

Corollary 1.3. Provided N ^> 1 and M 3> 1, in equation we may take 

(6) A = 0(iV- 1/2 + ikT 1 / 2 ) 

with an implicit constant dependent only on a and the implicit constants in the lower 
bounds on N and M. 



4 This paper contains a small error in line (0.4) that helped to spark the author's investigation 
into the van der Corput transform. 



5 A typo in the book has the logarithms in opposite places. 
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Corollary 1.3 allows us to take v as large as we like (and hence our intervals as 
large as we like) without increasing the bound on A. This, as well as the loss of the 
logarithmic factor, are common to applications of the results in this paper. 

More powerful results are possible. In many cases, the main theorems of this 
paper also allow one to extract the next term in the asymptotic for the van der 
Corput transform. Under the same general hypotheses of Theorems 1.1 and 1.2 we 
can improve the error terms to the following. 

Corollary 1.4. Suppose that f(x) and g(x) are real-valued functions with f four 
times continuously differentiable and g two times continuously differentiable on the 
interval [a, b}. Suppose there exist constants M , T , and U satisfying M = b — a 1, 
T ^> 1, and the bounds on lines ^ and 
Then, 

v^* / \ i fi \\ v^* 9(x r )e(f(x r ) -rx r + 1/8) 
2^ 9{n)e{f{n)) = ^ 



a<n<b f'(a)<r<f'(b) 



T(6) + T(a) + 0(-|=(l + f 



where T(/i) equals 

( g( rif {3) (rie(f(ri) _ g'(ti)e(f(n)) 

67ri/"(/x) 2 2mf"(/j) 
UM\ 



II/'(aOII=0 



°(^f) < ||r(M)ll < v^FCM) 



u v Tt i \\m\\>vm 



o 



M\\f>(x)\\^ M*\\f>(xW 



Here, [x] represents the nearest intege^to x; (x), the difference between x and the 
nearest integer to x, namely x — [x]; and \\x\\, the distance between x and the nearest 
integer to x, so that \\x\\ = \ (x)\. The function ip(x,e) equals 

1 .. e(rx) | . 1 

hm > tor e < 

R_>nr, Z-^i r _|_ f: 



2ni R^oo ' r + e 2 

0<|r|<_R 

which converges and is uniformly bounded for all real x. 

The size of the implicit constant in the big-0 term depends on the implicit constants 
in the relations of M, T, U , and the derivatives of f(x) and g(x). 



6 If x = n + 1/2 for an integer n, then it doesn't matter if we let [x] equal n or n + 1 provided 
that we do so consistently. 
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This in particular suggests that the size of T(/x) in Karatsuba and Korolev's result 

is optimal when | !/'(/•*) 1 1 > V 7"" (/■*)■ When ||/'(/x)|| = — that is, when /'(/i) * s an 
integer — the term T{jj) can bemore simply bounded by 0(UM/T). 



Corollary |1.4| can, in certain cases, be used to improve the Kusmin-Landau inequal- 
ity, a common result in the study of exponential sums. 

Theorem 1.5. (The Kusmin-Landau inequality — Theorem 2.1 in [9]^ 

Suppose f(x) is continuously differentiable and that f'(x) is monotonic on an in- 
terval [a, b\. Moreover suppose \\f'{x)\\ > 9 > on [a, b], where \\f'(x)\\ is the distance 
from f'(x) to the nearest integer. Then 

J2 e(/(n)) < cot (-° 

a<n<b ^ 

Corollary 1.6. Suppose f , T , and M satisfy the conditions of Corollary \l-4\ 

If 9 = min, eM ||/'(z)|| is positive, \\f\a)\\ > ^{a), and \\f{b)\\ > y/pf®, then 



e(/(6)) e(/(a)) 



2m(f'(b)) 2m(f'(a)) 

f 1 T If M 

+ V MO* + M 2 ® 3 + Vr Z \ 1 + 1 r 
where (x) represents the difference between x and the nearest integer to x. 



Corollary 1.6 follows immediately as a special case of Corollary 1.4| Moreover 



it suggests that the constant in the classical Kusmin-Landau inequality may not be 
optimal. The function cot(7r#/2) grows like 2/ir8 as 9 approaches zero; however, 



Corollary 1.6 suggests that the growth should be at worst like l/%9 as 9 approaches 
zero. 



While Corollaries L4 and L6 strengthen many previous results, they are con- 
strained to apply to short intervals. The results of this paper do not give similarly 
simple conditions and error terms when the size of the interval is allowed to grow 
arbitrarily large. But, given a specific sum, we can show very great improvements as 
one endpoint of the interval tends towards infinity. Consider the specific transform 

( 7 ) E* e ((if 2 ) = E* ( 24 -) 1/2 • e (- 4 - 3 + V8) + A, 

l<n<iV ^ ' (l/12) 1 / 2 <r<(V/12) 1 /2 

where N is an integer. (We have kept the (1/12) 1 / 2 due to it's natural appearance in 
the application of the van der Corput transform. It may be replaced by 1/2 with no 
change in value on the right-hand side, however.) 



If we apply the form of the error term from Theorem |1.1[ then we obtain the 
following result, which was included as an example in [23] . 

7 A short history of this theorem is given on page 20 of [S] ■ 
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Corollary 1.7. In line (J7| ; we have 

A 



0{N 



1/4N 



If we apply Theorem L2 from Karatsuba and Korolev instead, we obtain much 
finer results. 



Corollary 1.8. In line (J7| ; we have 



A 



0((logiV) 2 ) + I irmi 
0((logiV) 2 ) 



A/i/4 



\\{N/12)V*\ 



||(iV/12) 1 /2|| 
||(iV/12) 1 /2|| =0 



where ||x|| i/je distance from x to the nearest integer. 

In this case, the (logiV) 2 term comes from needing to break the sum on the left- 
hand side of ([7]) into roughly diadic intervals, each of which contributes an error term 
of size around logiV. One must also choose the endpoints of these intervals to be 
values which are 12 times a square, in order to keep the error term T(n) zero. 

Using the results of this paper, we can show that the result of Karatsuba and 
Korolev is almost sharp, in that we will extract an explicit term of size N 1/4 when 
iV 1 / 4 is smaller than || (iV/12) 1 / 2 ||~ 1 and extract an explicit term of size ||(A^/12) 1//2 || _1 
when ||(A^/12) 1 / 2 ||~ 1 is smaller than iV 1 / 4 . 

Corollary 1.9. In line we have 

' c + 0{N- l l 2 ) 
2?/;( v /iv7l^)(3iV) 1 / 4 e((iV/3) 3/2 + 1/8) 
+ 0(iV 3 / 20 ) 

+ 0(iV 5 / 12 ||(iV/12) 1 / 2 || 2 / 3 ) 



A 



(N/3) 



3/2 



x 



X 



27Ti((iV/12) 1 /2) 

+ c + 



||(iV/12) 1/2 || =0 

< \\{N/12)^ 2 \\ 
< (12iV)- 1/4 



1>(N, ((iV/12) 1 / 2 )) ||(iV/i2)V2|| > (UN)' 1 / 4 



m/^KN/uy/^w 

where c is a particular constan^ and ip(x) = ip(x,0) is the sawtooth function. The 
remaining functions are all as in Corollary l.J\ 

One way to understand these new error terms is via the geometry of the associated 
sums. Given functions / and g, one typically considers the curve S(t) : [0, 00) — > M. 2 
given by 

S(t) = J2 9(n)e(f(n)) + {t}g([t\ + l)e(f([t\ + 1)), 



Kn<t 



'Rough Mathcmatica calculations with N = 120, 000 give c 0.168 - 0.320i. 
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Figure 1. The curves associated to 
J2n<N e(nlogn) and it slightly extended 
1 vV 1 " 1 • e(- 



sn<N 

form ^ 



n<log(A r )- 



the exponential sum 
van der Corput trans- 

1 + 1/8) with N = 4000. 



where |_^J represents the floor of a real number x and {x} represents the fractional 
part of x. Geometric aspects of these curves have been well-studied [H El El [21] ■ 

Especially when /" is small, these curves generate a series of spiral-like figures. 
The van der Corput transform of the sum associated to S(t) then generates a new 
curve that can be seen to connect the center point of successive spirals by straight 
lines like in Figure [T} Thus the van der Corput transform can be seen as smoothing 
(if /" is small) or roughening (if /" is big) the curve. An easily accessible explanation 
for why this occurs is given in [I]. 

For our example case ^ , the corresponding curve is displayed in Figure [2j 

The term 



(8) 2^(v / A 7 7l2)(3^) 1/4 e((A^/3) 3/2 + l/8) 

from Corollary L9] describes the outer arm of the spirals and the connections between 
them. Figure [3 depicts the outer arm of the 500th spiral and its approximation using 
Corollary |1.9| The reason why this approximation appears so poor is due, in part, 
to the number of points on this part of the spiral. The mth spiral has on the order 
of m points in it. The outer arm of the spiral, where 

||(iV/i2) 1/2 || < (\2N)- l/ \ 

accounts for only about y/m points. The big-0 term O (iV 5//12 || (iV/i2) 1//2 || 2 / 3 ) is on 
the order of the explicit term g when || (TV/ 12) 1 / 2 1| x (12iV)~ 1/4 , so that g is a 
good estimate for only o(y/m) points, a negligible piece of the spiral. However, from 
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Figure 2. The curve associated to the exponential sums on the left- 
hand side of with N = 1,2, ... , 1200, and the 500th spiral of this 
curve — this and subsequent pictures are off-set closer to the origin. 



the picture, we can see that the approximation does a very good job of estimating 
the rotation around the spiral and appears to be of the correct order of magnitude 
for the distance to the center of the spiral. Although we do not calculate the big-0 
constants explicitly, it is quite possible the big-0 terms are smaller in modulus than 
(|8|) for ||(A^/12) 1 /2|| < (i2AT)" 1 /4_ 

It is possible to obtain far more explicit results along the outer spiral by making 
use of the Fresnel integral functions; however, this method does not provide good 
asymptotic data since it uses values of the Fresnel integral functions for which good 
asymptotics do not currently exist. More details will be given in Section [TTj 

The remaining points (in fact, almost all of the points) are described by the term 

< 9 > e 'W 3 ) 3/2 )( 2 ri (Wl2)^) -^<W 12 ' 1/2 »)- 

which gives the tight inner weave of a given spiral. The term ip(N, ((iV/12) 1 / 2 )) can 
be bounded by 0(1) for all N. The big-0 term 0{N- l ' 2 \\{N /12) 1 ' 2 W^) is bigger 
than 0(1) for only 0(m 2 ^ 3 ) of the points on the mth spiral, again a negligible piece. 
Thus, we know almost all points on the inner weave (and hence almost all points on 
the full curve) to within a 0(1) error and often much better. 

The big-0 terms grow to the size of the explicit terms Q and @ as \\{N /Uf^W 
approaches (12A^)~ 1 / 4 in both cases, so we don't have good information about how 
the curve transitions from the inner weave to the outer arm and vice-versa. We also 
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Figure 3. The outer arm of the 500th spiral of the curve associated 
to §f§ and its approximation by Corollary 1.9 





Figure 4. The inner weave of the 500th spiral of the curve associated 
to ([7]) and its approximation by Corollary 1.9 ignoring the contribution 
of the ip(x, e) term. 



do not have a good heuristic explanation for why e(l/8) appears in (|8|) on the outer 
arm, but \ji = e(— 1/4) appears in ^ for the inner weave. 
This paper is arranged as follows. 
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In Section [2j we will provide a short catalog of notations that are frequently used 
in this paper. 

In Section |3j we will outline the methods of this paper and how they overcome 
difficulties encountered in other results. 

In Section |4j we will state the main theorems of this paper as well as a briefly discuss 



them and their assumptions; in this section, we also include a proof of Corollaries |1.3 
andO 

In Section [5j we will cite a number of lemmas from other papers and books that 
will be necessary in the proof of the theorems, as well as give propositions that show 
how the lemmas apply under the specific conditions of the main theorems. 

In the remaining sections, we prove the primary theorems of our paper and then 
prove Corollary |1.9| 

2. Notation 

We will frequently use the Landau and Vinogradov asymptotic notations. The 
big-0 notation f(x) = 0(g(x)) (equivalently, f(x) <C g(x)) means that there exists 
some constant c such that |/(x)| < c|g(:r)| on the domain in question. By 0(f(x)) = 
0(g(x)), we mean that a function which is asymptotically bounded by f(x) will 
also be asymptotically bounded by g(x). The little-o notation f(x) = o(g(x)) as 
x — > a — typically with a equal to oo — means that 

lira 44 = 0- 

g(x) 

By f(x) x g(x), we shall mean that g(x) <C f(x) <C g{x). 

We will need several functions that relate a real number x to the integers nearby 
it. Let [x\ denote the usual floor of a real number x, the largest integer less than 
or equal to x, and let \x] denote the usual ceiling of a real number x, the smallest 
integer greater than or equal to x. Let {x} := x — [x\ denote the fractional part of 
a real number x. Let s(x) := {x} — 1/2 denote the sawtooth function, and let 



ip(x) 



s(x) x $l Z 

o x e z 

denote the "smoothed" sawtooth function, with Fourier series 



OO 



, 1 sin(27rra;) 
W = — V • 

r=l 

We also wish to have a modified sawtooth function, given by 

1 e(rx) | 1 

ib(x,€) := hm > for e < -. 

0<\r\<R 



The convergence of ip(x, e) will be guaranteed by Proposition 5.16 
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Figure 5. The real (left) and imaginary (right) parts of the function 



Let (x) := s(x + 1/2) denote the difference between x and the nearest integer to x, 
let [x] := x — (x) denote the nearest integer to x, and let \\x\\ := min{l — {x}, {x}} 
denote the distance from x to [x]. Let be given by 



x 



\x\ 



\x\\ ^ 
\x\\ = 



A few more definitions will simplify the (nonetheless still complicated) statement 
of the theorem. For visibility's sake, we will avoid writing (x) all the time when the 
choice of argument is always the same. Given functions / and g, let 

H = gf {3) + 3g'f" 

G = \2gg'\ff 

(2g"Yg> (2g") 3 (f"g) 



W. 



± 



Wn 



r± = f 
r = f 



(H ± VH 2 - G) 2 (H ± VH 2 - G) 3 
H 2 f^ 

27g(f") 5 

, H ± VH 2 - G 



2g" 

, Uf"? 
H ' 
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and given an integer r, let 

(/ # (x)-r)^(x)- g (x)r(x) 

3. Heuristics for the van der Corput transform and the new ideas 

of this paper 

The transform is a simple two-step process: the Poisson summation formula is 
first applied to obtain a sum of integrals, and then each integral is estimated using 
stationary phase methods. In particular, the Poisson summation formula (see Lemma 



5.15) gives 

* R r b 

(10) S~] g(n)e(f(n)) = lim S~] / g(x)e(f(x) - rx) dx. 

a<n<b r =-R Ja 

An oscillatory integral, like 

(11) I g(x)e(f(x)) dx, 

J a 

is said to have a stationary phase point if there exists an x' G [a, b] such that f'(x') = 
0; if such a point exists, then we expect (11) to be roughly 

g(x>)e(f(x>) + 1/8) 

plus some small error, where c(x') equals 1/2 if x' equals a or b and 1 otherwise. In 
our case, the stationary phase point of the integrand 

g(x)e(f(x) - rx) 

occurs at x = x r , which is defined by f'{x r ) — r = 0. Provided f'(a) < r < /'(&), this 
stationary phase point will be inside the interval of integration, and so we expect 

v^* / w ,/ u v^* 9(x r )e(f(x r ) - rx r + 1/8) 
2^ 9(n)e{f(n)) w ^ T^rl ' 

a<n<b f'(a)<r<f'(b) V/ \ Xr > 

The error term for the van der Corput transform thus comes from the error in 
estimating all the integrals of (10). 

The first problem in doing such an estimation comes from the integrals without a 
stationary phase point. Typically the method one uses to estimate such an integral 
is an application of integration by parts. 



^ c(f(x) — Tx\ 

e(f(x) — rx) dx — 



f"( x ) 

+ / o 7777 S \2 e (f( x ) ~ rx ) dx 

a 27Tl(/'(x) - r) 2 



27ri(f'(x) - r) 

(We will assume that g(x) = 1 for a while to simplify the heuristics.) The terms 
e(f(x) — rx)/2iri(f'(x) — r) at a and b are referred to as the first-order endpoint 
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contributions and exhibit a great deal of cancellation (see Proposition 5.16); for 
example, if f'{a) is an integer, then 

and similarly if /'(&) is an integer. 

However, we are still left with the integral 

f"(x) 



-e(/(x) -rx) cte 
2m{j (x) — r) z 

to evaluate. We could estimate the integral by taking absolute values of the integrand, 
but we would then obtain an estimate of the size 0((r — /'(a)) -1 ) +0((f'(b) — r) _1 ), 
which is of the same order as the first-order endpoint contributions, but without the 
cancellation that would allow them to be nicely summed. 

One possible solution would be to apply integration by parts a second time: this 
gives 

(12) f\ jy} .ifc)-,)*-^'-"' 11 



2m(f' (x) - r) 2 ww ' (27ri) 2 (f (x) - r) 3 



e(f(x) — rx) dx, 



dx \(27ri) 2 (f'(x) -r)\ 

with a second-order endpoint contribution and a new integral. The second-order 
endpoint terms are absolutely convergent, and if f(a) and f'(b) are integers, they 
sum to 0(f"(a)) + 0(f"(b)), which is a good estimate if /" is small at the endpoints. 



On the other hand, if we try to bound the new integral in (12) by taking absolute 
values of the integrand, we get terms of the same order of magnitude as the second- 
order endpoint contributions (which, unlike the first-order endpoints, are summable). 
But in addition, we have the total variation of f"(x)/(f'(x) — r) 3 on [a, b], which 
may be roughly bounded by the sum of the moduli of local maxima and minima of 
f"(x)/(f(x)-rf on [a,b\. 

Unless f(x) is quadratic — as in Coutsias and Kazarinoff's case^j — it is difficult to 
find good estimates on these integrals. More applications of integration by parts or 
the presence of a non-constant g only make things worse. Therefore, these terms are 
often avoided entirely by applying the truncated Poisson formula instead of the full 
Poisson summation formula. 

Proposition 3.1. (Truncated Poisson formula — Proposition 8.7 in [13\) p°| 

9 Coutsias and Kazarinoff actually analyze the resulting integrals and endpoint contributions in 
the quadratic case when inegration by parts is repeated many times. 

10 An explicit version of the truncated Poisson formula for non-trivial g is given in Lemma 7 of 
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Let f(x) be a real function with f"(x) > on the interval [a,b]. We then have 
(13) e (/M) = E [ <f( x ) - rx ) dx + °( e_1 + M/? - « + 2)) 

a<n<b a-e<r</3+€ a 

where a, (3, and e are any numbers with a < f'(a) < f'{b) < /3 and < e < 1, i/ie 
implied constant being absolute. 

The truncated Poisson formula is the source of the 0(log(/'(6) — f'(a) + 2)) error 
term in many of the results mentioned above. 

After applying the truncated Poisson formula with /3 = f'(b) and a = f'{a), one 
is left with approximately f'(b) — f'{a) integrals with which one hopes to apply 
stationary phase estimates. These estimates work best on a small interval around 
the stationary phase point, where the second derivative of / is fairly constant, and 



the higher derivatives of / are small. So we would like to break the integrals in (13) 
into several pieces, such as 

(14) 

where the point of stationary phase is near the middle of [a 1 , b'\. However, to make 
this effective, we would again require good estimates on integrals with no stationary 
phase point (the integrals from a to a 1 and b' to b). In addition, one is faced with 
possible first-order (and higher) endpoint contributions at a' and b' . (While estimates 
of stationary phase integrals benefit from integrating on a small interval, they suffer 
again if the interval is too small. The 0{T(a) +T(6)) terms on line Q that appear in 
stronger results on the van der Corput transform arise from bounding the stationary 
phase integrals where x r is closest to a or 6.) 

Because of this, many results on the van der Corput transform seek to treat the 
entire integral from a to 6 as a single stationary phase integral, hence the common 
conditions, as in Theorem |1.1| that f"(x) has constant order on the entire interval 
[a, b] and that higher derivatives of / be likewise small on the entire interval [a, b]. 

The techniques of this paper seek to overcome some of these difficulties. 

First, we use the recent and powerful stationary phase estimates first proved by 
Huxley [HI H2] and refined by Redouaby and Sargos [30] . These estimates show that 
stationary phase integrals, such as the middle term in ( Jl4| ), also generate first-order 
endpoint contributions. In fact, these contributions directly cancel the first-order 



endpoint contributions at a' and b' generated by the first and third term in (14). 
Thus we no longer need to treat the full integral from a to b with stationary phase 
estimates and so can replace the global restrictions on the derivatives of / and g 



(such as those in Theorem 1.2) with local restrictions. 



In particular, we replace the the standard assumptions 

f"(x) x TM~\ 
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\f (3) (x)\ <TM~ 3 , and 
|/< 4 >(i)|«TM- 4 
for x G [a, b], with an assumption that looks like 

f"(z) x T{x)M{x)-\ 
\f( 3 \z)\ <£T(x)M(x)- 3 , and 
|/ (4) (^)| <T(a;)M(a;)- 4 



for z G [x — M(x),x + M(x)\ and x G [a, 6] with some function Mix) (see Section 4.1 
for more details about the function M(x)). 

Second, to avoid use of the truncated Poisson summation formula, we develop a 
method to get reasonable bounds on the integrals arising from applying integration by 



parts twice, as in (12). Instead of looking at one integer r at a time and counting the 
contribution of f"(x)/(f'(x) — r) 3 for each x's which give local maxima and minima, 
we instead look at each x and imagine a possible fixed real-valued r that causes the 
point x to be a critical point of this function. This then defines a function r(x), and 
we estimate the sum over all values of x where the function r(x) takes integer values 



using a variant of Euler-Maclaurin summation (see Proposition 5.14 and Section 6.5). 

Third, since the previous two techniques frequently benefit from /" being large, we 
develop a method to deal with the large second-order endpoint contributions in this 
case. The reason for these terms being large is that 

/» 
(/'(a) - r)3 

is large when r is close to f'(a) (and likewise at b). However, if f"(a) is large and 
f^(a) not too large in comparison, then a small shift in a to, say, a + e should greatly 
increase the size of the denominator while keeping the numerator roughly the same 
size. Therefore, by altering the endpoints of the integrals in the Poisson summation 
formula by a small amount when r is close to f'(a) or /'(&), we can reduce the size 
of the second-order endpoint contributions. 

With these techniques, the main theorems of this paper give the van der Corput 
transform in the following form. 

v^* / \ ( £ ( \\ ^* g{x r )e{f{x r )-rx r + \) 
2^ 9{n)e{f{n)) = ^ -=== «- 

a<n<b f'(a)<r<f'(b) V J ^ %r ' 

- V{b) + V(a) + 0(A 1 + A 2 + A 3 + A 4 ) 

The T> terms will usually equal the first-order endpoint contributions and will be 
explicit in many cases. The term Ai will estimate the error that arises when there 
exists an x r close to a or b and is analogous to the T(a) + T(b) terms present in the 
results of Kolesnik, Liu, and Karatsuba and Korolev. The term A 2 will estimate the 
second-order endpoint contributions at a and b. The term A 3 (and a bit of A 4 ) will 
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estimate the remaining error in applying stationary-phase methods and the second- 
order endpoint contributions at a 1 and b'. The remainder of A4 will estimate the size 
of the integrals that arise from applying integration by parts twice. 

The techniques outlined in this section allow for a number of interesting extensions 
to the van der Corput transform. 

Since we replace the global restrictions on the derivatives of / and g with local 
restrictions, we can apply the van der Corput transform to many new sums, including 
those where / and g have moderate oscillations. Likewise, we can now directly apply 
the van der Corput transform to sums where /" is large. 

The function f{f'~ 1 {x)) — xf'~ 1 {x) that appears on the right-hand side of the van 
der Corput transform is sometimes referred to as the van der Corput reciprocal 
of the function /. If the second derivative of / is large, then the second deriva- 
tive of the van der Corput reciprocal of / will be small, so the results cited above, 
if they cannot be applied directly to a sum, can be — and frequently are — applied 
"backwards" to the van der Corput transform of the sum. 



For example, if we attempted to apply the form of the error from Theorem 1.2 
directly to the sum 

0<n<N 

for some (3 > 1, then the error term 0(log(/'(6) — f'(a) + 2)) would be at least O(N), 
the size of the trivial bound on the sum (and that doesn't even take into account the 



size of the implicit constant!). Theorem 4.1 allows us to apply the transform to this 
sum directly and obtain a 0(1) error. 

In general, if /" > 1 is reasonably large compared to g and both are free of wild 
oscillations, then A should be no larger than 0(max a < x <6 g(x)). 

As a final note, the results of this paper do not give an improvement in all possible 
cases. If, for example, /" varies between a very large value and a value very close 
to zero, then our results may give very poor error terms. Another case where our 
results do not give any improvement is 

sm(an)e(/3n 2 ). 

a<n<b 

4. Details of the theorem 

4.1. The auxiliary function M(x) and condition (M). We will use a function 
M(x) to measure the length of an interval around x where the functions f"{x) and 
g(x) are well-approximated by their Taylor polynomials up to the second degree. The 
larger M(x) is, the more linear the functions f"(x) and g(x) will appear at the point 
x. 



11 Redouaby and Sargos [31] and the papers cited within have studied the van der Corput recip- 
rocal in much greater detail, including useful asymptotics for how the reciprocal changes as / is 
perturbed. 
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To be more concrete, by condition (M), we shall refer to the existence of positive, 
bounded, continuously differentiable functions M(x) and U(x) on [a, b], along with 
several associated positive constants C 2 , C 2 -, C4, D , Di, D 2 , 5, which together 
satisfy several conditions: 

(I) max{M(a),M(6)} < b - a; 
(II) 5 < 1, and let r] := 35/C 2 - satisfy r] < 2; 
(III) If J:=[a- c{a) ■ M(a),b + c(b) ■ M(b)}, where 



c(x) 



0, if Z n (f'(x) - f"(x), f(x) + /"(*)) \ {/'(x)} = 0, and 

1, otherwise, 



then / is C A [J] and g is C 3 [J]; and, 
(IV) If I x denotes the intersection [x — M(x), x + M(x)] D J, then for all x G [a, b] 
we expect all the following inequalities to hold for z G I x : 

1 -f'(x) < f"{z) < C 2 f'(x), \g(z)\ < D U(x), 



2- 



l/ (4) WI<^4^4, and |/(^| <D 2 U{,) 



M(x) 2 ' ISW| - z M(a;) 2 ' 

For any given /(x), ^(x), a, and 6, there are infinitely many possible choices of the 
functions M(x), U(x), and the associated constants that satisfy the above conditions. 
Therefore, for the remainder of the paper, we shall assume that if f(x), g(x), a, 
and b remain unchanged, the particular functions M(x) and U (x) and associated 
constants we reference will be likewise unchanged. If f(x) and g(x) are unchanged, 
but a and b allowed to vary, then we will assume, first, that the associated constants 
will be unchanged, and, second, that given two intervals [ai,foi] and [02,62] with 
corresponding auxiliary functions M±(x), U±(x) and M 2 (x), U 2 (x), we have M\(x) = 
M 2 (x) and Ui(x) = U 2 (x) on (a±,bi) fl (a 2 ,b 2 ). 

Since condition (M) is rather intricate, we pause a moment to illuminate it further. 

(II) The requirement of 5 < 1 and r] < 2 in condition (M) part (II) guarantee that 



the Taylor approximations to /" has good properties (see Proposition 5.10 



Lemma 5.1, and the remark following the lemma). 
(Ill) The restriction of condition (M) part (III) to have / and g be several times 
continuously differentiable on the larger interval J is often no restriction at 
all. Many applications of the van der Corput transform have f(x) and g(x) 
be polynomials, exponentials, logarithms, or other C°° functions that exist on 
large domains. 

The bounds on the derivatives of /" and g in condition (M) part (IV) allow 
us to use estimates on stationary phase integrals that will be discussed in 
Section [5j By expecting the properties of /" and g to extend to the larger 
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interval J, we can also extend the integrals under consideration to this larger 
interval. 

(While typical stationary phase results require only that g is twice continu- 
ously differentiable, we require three times in order for the functions W'± and 



r± to exist.) 

(IV) To understand the complex system of inequalities in condition (M) part (IV), 
it is helpful to rewrite them in terms of Taylor approximations. The first three 
inequalities then become 

f"(z) x f"(x) 

^ = ^( 1+ °(#- a) 

(15) f"(z) = f"(x) (l + O (^-Az - xf ) ) + f&>(x)(z - x) 



M(x) 

for z G I x with implicit constant 1 in the second and third lines. The equality 



on line (15) gives definite form to our earlier statement that the larger M{x) 
can be, the more linear /" appears at the point x. 

Alternately, one can interpret M(x) as somehow representing the rate of 
decay as one takes successive derivatives. At the point z — x, the various 
inequalities imply — among other things — that |/( 3 )(a;)| <C f"(x)/M(x) and 
|/( 4 )(*)|«f(x)/M(x) 2 . 

The conspicuous absense of a constant C 3 is intentional. We could include 
such a constant, but later estimates are made slightly simpler by having C3 = 
1 

Here are some examples of the function M(x) in different circumstances. In each 
case we assume that the associated constants always take the same values: 6 equals 
1/2, and C2, C 2 -, C4, D , D\, and D 2 all equal 2. 

• If f(x) and g(x) are polynomials, there exists an e > such that condition 
(M) is satisfied with M(x) = e ■ x and U(x) = g(x) for sufficiently large x. 

• If f(x) and g(x) are exponential functions (a- (3 X , with (3 > 1), there exists an 
e > such that condition (M) is satisfied with M(x) = e and U(x) = g(x). 

• If f(x) = t(\ogx)/2ir and g(x) = x~ a , which corresponds to the Riemann 
Zeta function at z = a + it, then there exists an e > such that condition 
(M) is satisfied with M(x) = ex and U(x) = g(x) independent of t. 

• If f(x) = ax 2 + sin(7a;) and g(x) = 1, then there exists an e > such 
that condition (M) is satisfied with M(x) = e^fx and U(x) = 1. 

• If f(x) is a power function and g(x) = sin(ax) for some constant a, then 
there exists an e > such that condition (M) is satisfied with M(x) = e and 
U(x) = 1. 
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4.2. The assumptions. Assume that condition (M) holds for some function M(x) 
and U(x). Let J be as in condition (M) part (III). 

We assume that f"(x) is a positive, real- valued function, bounded away from 
on the interval J, and that g(x) is a real-valued function on J as well. (The case 
f"(x) < may be considered by taking the conjugate of the sum.) This will guarantee 
that f'(x) is a continuous, monotonic function, so that x r := / /-1 (r) is well defined. 

Let J± be the union of all intervals [a', b'] C J such that G(x) ^ and H(x) 2 — 
G{x) > for x G [a', b'], let Jo be the union of all intervals [a',b'] C J such that 
g"(x) = 0, g(x) ^ 0, and H(x) ^ for x 6 K,&']? and let J be the set of points 
x G J such that g(x) = but g'(x) ^ and g"{x) ^ 0. It is possible that J± will 
contain isolated points due to H(x) 2 — G(x) having a on an interval where it is 
otherwise negative, and it is also possible that J will contain isoluted points due 
to g"(x) having a on an interval where it is otherwise non-zero; we denote the set 
of isolated points in J± and J by J± and Jq respectively. Let dJ± and dJ denote 
the endpoints of the non-zero-length intervals contained in their respective sets. In 
particular, dJ± nJ| = and dJ Pi Jq = 0. 

For our final assumption, suppose that if (a', b') is an interval contained in J±, then 
either H(x) 2 — G(x) equals on the whole interval and H(x) does not equal at any 
point on the interval nor does it tend to at the endpoints, or H(x) 2 — G(x) does 
not equal at any point on the interval nor does it tend to at the endpoints; and 
if (a', b') is an interval contained in Jo, then g(x) does not tend to at the endpoints 
of this interval. 



4.3. Statement of the theorems and proofs of some corollaries. 

Theorem 4.1. Assume all the conditions of Section \4-£\ hold. Then 

* — ■ \* \* q(x r )e( f (x r ) — rx r + I) 

9{n)e{f{n)) = £ 9[ r) ^=g= ^ - 7>{b) + V(a) + A 



a<n<b f'(a)<r<f'(b) 

where 

V(x) = V°(x)+V*(x) 

( g(x)e(f(x) + lf{x)jx)x 



V°(x) 



( 1 , ,t I ft ^\^^ f"( x ) < Wf'( x )W 

g(x)e(f(x) + lf(x)j x)^(x, (f'(x))) \\f'(x)\\ < f"(x) < 1 - 

0(U{x)) l-\\f'(x)\\<f"(x)<l 

O I U(X) I 1 + — ^- + ; 1 I I fix) > 1 
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g (x)f^(x)e(f(x)) g'(x)e(f(x)) ||f , M||=n 

Qmf"(x) 2 2mf"(x) U 1 J " 

.0 \\f'(x)\\^0 



^0(A i (a) + A,(6)) + A 4 



i=i 



^4=,^t^!> ll/'Mll^o. ,,,,.>] 



C/(x) 



\\f'(x)\\=0 



f"(x) 2 (b - a) 3 

\\f(x)\\^0, mx = 

f U(x, . , 
mm < 1 



+ < 



M{X) 1 I{X) > i 1 i/||/'(x)||=0 O rm,> 
+ ^)mm(r(,),_) 

C/(x) U(x)f"(x) 

+ „ xL otherwise 



{M(x)\\f'(xW \\f'(xW 
? ' ^ V ' ^ ■ ' ,) <fs 



s - a) 3 V /"(x)M(x) - a) 

+ f»(a)*(a-a)* + f»(b)*(b-a) 
U(x) 

a f"{x){b-xf y ' f"(x)M(x) ' /"(x)(6-a;) 

U(b) U(a) 
+ f"(ly(b-b) 3 + f"{af{b-af 



i + v^5)m(x)) (1 + -; 



/' '(:r)M(:r) 3 V v J v ' v 7 ^ /"(x)M(x) 
+ /C( Jo, W , r ) + /C( J±, W + , r+) + /C( J±, r_) 

^ g'(x)f"(x) 2 

J null 



dx 



JC(I,W,r)= (\W(x)\\r'(x)\ + \W'(x)\) dx+^\W{x) 
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+ \s(x)-W(x)\. 

xGl and r' changes sign at x 
or x£dl 

The number m x equals the cardinality of the set 

Zn(f(x)-f(x)J'(x)+f"(x))\{f(x)}. 

In particular, 

'0 f"(x)<\\f'(x)\\ 
0(1 + /"(*)) f"{x) > \\f{x)\\ 

Also, a equals the smallest value in the interval [a + min{M(a), C^ 1 }, b] such that 
/'(a) is an integer; if no such value exists, then we may take A 3 (a) = 0. Similarly, b 
equals the largest value in the interval [a,b — min{M(6), C^ 1 }] such that f'(b) is an 
integer, and if no such value exists, then we may take A 3 (5) = 0. 

The functions M(x) and U(x) are as in condition (M) ; the sets J±, J±, Jo, and 



Jo are all as in Section 4-2, and the functions ip, s, Wo, r , W±, r±, (■), [■], and || • || 
are all as in Section\E 

The implicit constants in the big-0 terms are dependent upon the constants C2, 
C2- , C4, Do, D\, D2, 5 in condition (M). The convergence of the integrals in A4 



will be guaranteed by Proposition 5.17\ 



Before continuing to further refinements of the main theorem, let us pause to 
assure the reader that the vast collection of error terms above are more psychologically 
daunting than they are computationally difficult. As an instructive example, consider 



the case of Corollary 1.3 where 



'<•>-!(£)■. 

a = N, and b = vN . N and X are assumed to be positive and a and v both strictly 
greater than 1. Without loss of generality, we may assume N < \/X (otherwise, we 
would instead consider the right-hand side of ((5j)). 

Note that f(x) and g(x) are both power functions, so condition (M) holds with 
C2, C 2 -, C4, D , Di, and D 2 all equal to 2, 5 equal to 1/2, M(x) equal to ex for some 
small e > dependent on a only, and U(x) equal to g(x). 

The benefit of our conditions 1 N < \fX is that on the interval [N, uN] , the 
functions f"(x) and M[x) are both bounded below by some constant (depending 
on a and the implicit constant in 1 N). This immediately gives that the T>(x), 
Ai(x), and A 2 (x) terms are all 0{U{a) + U{b)) = 0(N^ 1 / 2 ), since each summand of 
these three terms contains a factor of U(x) divided by some positive power of f"(x) 
and M(x), which are bounded from below. (Note that ||/'(x)|| terms only appear if 
||/' (a?) || is bigger than f"(x) to begin with.) 

For the A 3 terms, we need to consider a and b. It might be quite difficult to 
calculate these values explicitly, but we only make the A 3 terms bigger by assuming 
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a is as small as possible (that is, a = a + min{M(a), }) and b is as large as 
possible. Since M(x) > 1 on [N, uN], we have that a — a ^> 1 and b — b^> 1. Thus, 
by first applying the fact that /"(a;) 3> 1 and M(x) ^> 1 and then that a — a ^> 1, 
the integral in A3 (a) may be bounded as 



U(x) 



1 1 

1+ ,,„ r~r~T7 ; + 



ff /"(x)(x - a)3 V /"(x)M(x) /" (x)(x - a) 



dx 



< [/(a) 



1 1 

+ — ax 



o (x-a) 3 (x-a) 4 



C/(a) 



a — a) 



< t/(a) 



which is 0(iV~ 1//2 ) again. A similar argument shows that the remaining terms in 
A 3 (a) and all the terms in A3 (6) are 0(N~ 1 ^ 2 ) as well. 

This is one reason why the error terms are much easier to deal with than at first 
glance: provided f"(x) ^> 1 and M(x) > 1 on [a, b\, the error terms T>(x) and Aj(x) 
for i = 1,2,3 are all bounded by 0(max[ a) b]U(x)), where the implicit constant is 
dependent on the implicit constants in the lower bounds of f" and M. 

To nicely bound the first integral in A4, however, we must use the fact that M(x) 
and f"(x) are not only bounded below, but also somewhat large. In particular, we 
have 



f"(x)M(x) 



3 / 2 dx 



1 + 



1 + |M'(x)| 
f"(x)M(x) 



dx 



la M(X 

poo 

< / X 
J N 

< N' 1/2 



The greatest complication comes in evaluating the remaining condition from Sec- 
tion ^2 dealing with the functions G and H and the remaining error terms in A4. 
In this case, the function g" is never zero on the interval J = [N(l — e), uN(l + e)]. 
Hence Jo and J nu ii are empty, so both )C(Jo,Wo,r ) and the final sum of A 4 are 
0. This also shows that the functions / and g satisfy the remaining condition from 
Section gJJ 

Since / and g are power functions, we can quickly verify, with simple hand calcu- 
lations, that the functions H, G, W + , and W_ are all power functions as well; that 
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is, each is of the form c\ ■ x C2 . In particular, the functions W±(x) are equal to 

where c aj ± is a constant depending on a and the sign of ±. The functions r±(x) are 
equal to 

C n -i-~rz X ~\~ C n _i_—-z X , 

with new constants dependent on a and the sign of ±. 

Also, H(x) 2 —G(x) is a power function; so, J± is either the interval [N(l— e), uN(l+ 
e)] or is empty (if H(x) 2 — G(x) < for all x, in which case the remaining /C terms 
equal and we are finished). 

In general, we expect that if g(x) and f"{x) share the same rate of decay in their 
derivatives, then the various functions H, G, and the W and r functions should all be 
relatively well-behaved. In our case, since g and /" are both power functions, taking 
derivatives of each entails multiplying by a constant and dividing by x. Thus, we 
have that g(x)f^(x) is a constant multiple of g'(x)f"(x), hence why if is a power 
function as well. 

Returning our attention to the remaining error terms of A 4 , first consider the 
integrals. We have 

/ {\W4x)\\r'_(x)\ + \W + (x)\\r' + (x)\) + {\W' + (x)\ + \W'_(x)\) dx 
Jj± 



nvN(l+e) jya i j^2a l 

< / — X~" a + x~2~ 2a dx 

JN(l-e) X X 

N l/2 

The functions r' + (x) and r'_(x) have at most 5 zeroes on [a, b], and since the functions 
W + (x) and W^(x) are decreasing and since J± is empty, the remaining terms in A 4 
are bounded by 

|s(iV(l - e)) • W-(N(1 - e))| + \s(N(l - e)) • W + (N(1 - e))| 



< N~ 1/2 . 



(We implicitly used that iV » 1 and M » 1 implies X > 1.) 
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Since we assumed N < \/X, we have N^ 1 / 2 > M -1 / 2 , this completes the proof of 



Corollary 1.3 



Under fairly strong and yet quite common conditions, we can remove much of the 
complication from the assumptions of Section <L2_ and from A4. 

Theorem 4.2. Assume that f and g be real-valued functions satisfying f"(x) > on 
[a, b] and condition (M) for a function M(x) such that 

M(x) > max{6 — x, x — a} 

on [a,b]. Also assume that m a = = 0, where, as in Theorem J^.l, m x is the 
cardinality of the set 

Zn(f'(x)-f(x),f(x)+f"(x))\{f(x)}. 



Under these conditions, the result of Theorem 4-1 holds, with A 4 now just equal to 



U(x) 



1 + y r rix)M(x / 



1 



1 + \M'{x) 
f"(x)M(x) 



dx. 



I a f"{x)M{xf 

The condition that m a = mb = also simplifies more of the terms to the following. 

{ g(x)e(f(x) + lf(x)}x)x 

x 



V(x) + 0{A 1 (x)) 



g(x)f^(x)e(f(x)) g'(x)e(f(x)) 



6nif"(x) 



+ 



U(x) 



2nif"(x) 



f"(x) 2 (b-ay 



\\f'(x)\\=0 



Using Theorem 4.2, we may now prove Corollary |1.4[ 

For Corollary 1.4| recall that the conditions are that f(x) and g(x) are real- valued 
functions, f(x) four times continuously differentiate le and g(x) two times continuously 
differentiable on [a, b\. Also, 

f"{x) > T/M 2 f (2+r \x) < T/M 2+r g^ r \x) < U/M r 

for r = 0, 1, 2 on [a, b], for constants T>1, M = 6- a>l, and U. 
First, we wish to remove from consideration the case when 



(is) o<n/v)ii< Vr¥) 

for fj, equal to a or b. 

Suppose that not only is (16) satisfied for both a and b, but that |/'( a )l = 
that is, the nearest integer to f'{a) and f'(b) are the same. In this case, we have by the 
mean value theorem that f'(b) — f'(a) = (b — a)/"(0 for some £ G [a, b}. But by our 
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assumption on the size of ||/'(a)|| an d we nave that fib) — f'(a) <C y/T /M. 

Combining these we have that b - a < f'ffl^y/T/M < M/y/T. Thus we have 

'UM\ 



J2*9(n)e(f(n)) = o( 1 ^) 
<r,<h V V J / 



a<n<b 



We can add in the sum 



E* 

f'(a)<r<f(b) 



g{x r 



--e{f{x r ) - rx r + 1/8) 



to the right-hand side, since this sum by assumption has at most one term, which is 
of size UM/y/T. Thus the corollary is proved in this case, since U/y/T ad UM/T 3 ^ 2 
are bounded by UM/ y/T. 

If [/'(a)] does not equal [/'(&)], then we may replace an endpoint /i where (16) 
holds with fj,' the closest value to // inside the interval [a, b] where ||/' (/•*') II = \J /"(A*') 
or || /'(//) || = 0; by the same argument as the previous paragraph, doing this generates 
an error of size 0{U M / y/T). The only way we could not do this is if [/'(a)] + 1 = 
[/'(&)] and f"(x) is at least 1/4 at some point x on [a, b], in which case, f'(b) — f'(a) < 
1 <C y/T /M, so that b — a <^ M/y/T and the corollary follows in this case by the 
same argument as in the last paragraph. 



Thus it suffices now to prove Corollary |L4 in the case where H/XaOII — \J /"(a-0 
or ||/'(aOII = f° r ^ equal to a and b. The conditions of Theorem 4.2 hold with 
M(x) = b — a on [a, b]. The size of the constants in condition (M) occur in the 
bounds on the derivatives of / and g in terms of the constants T, M, and U, and 
U (x) just equals U. 

Now, we go through the error terms of Theorems 4.1 and 4.2 First, we have that 

+ 1>(n, (fin))) 



£(//) + O(A x (a0) 



II/'CaOII > v/FR 



67Ti/"(/i) 2 



ll/'(*)ll=0 



+ 



V T 2 I 

For the A2(x) terms, we use that M(x 

UM 



A 2 (/i) < 



y2 
+ 



(i + Vt) (i 



M and f"{x) 
T \ 



T/M 2 to obtain 



r 



U 



M||/'(/.)|| 2 
UM U 



+ 



m 2 y 

UT 



M 2 \\m\\ 3 

UT 



T3/ 2 M||f(/i)|| 2 M 2 ||f(/i)||3- 
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For the A3 (a;) terms, we again need to understand a and b. In the proof of Corollary 



1.3, we used the fact that /" and M were generally quite large to get good bounds. 



Here instead, when /" is small, we shall show that a — a and b — b are going to be very 
large most of the time. In particular, by the mean value theorem, for some £ G [a, b] 
we have a-a = (/'(a)-/'(a))//"(f) > \\f(a)\\M 2 /T since /'(a) must be an integer. 
Thus for A 3 (a), we have the following bound. 

3{a) = j a r{x){x-af \ + f"(x)M(x) + f"(x)(x-a) ' " 
+ 9(a) + 9(b) 



f"(a) 2 (a-a) 3 f"(b) 2 (b-a) 3 
UM 2 ( M M 2 \ UM 4 UM 

< 7777= 777 1 + 777" + 7777= 7 + 77777= 77 + 



T(a-a) 2 \ T T(a- a) J T 2 {a - af T 2 
UT / M 1 \ C/M 



M2||f(a)P V T ||/'(a)||; T* 



<<C M||/'(a)P + M 2 \\f'(aW + ~T 2 ~ 

A similar bound holds for A3 (6). 

Now for the last term A 4 , we have the following much simpler calculation. 

Ja f"{x)M{xf \ ' V f"(x)M(x)J 



U ( M 

« -F= 1 + 



/TV T 
This completes the proof of Corollary |1.4 



We conclude with two additional theorems which extend the results of Theorem 



4.1 in specific directions. The following theorem shows that if the various error terms 



from Theorem |4.1| are nicely bounded as b tends to infinity, then many of them can 
be replaced by a constant plus a o(l) error. 



Theorem 4.3. Assume all the conditions of Section ^2 hold with a fixed and b 
tending to infinity. 

Let Kb be the set of all x G [a, b] such that x + M(x) > b, where M(x) is as in 
condition (M). Let dS for a set S have the usual meaning of the boundary of the set 
(the endpoints of all intervals and the isolated points). 

Suppose we have the following functions 
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U(x) 



K h 



1 + 



+ 



nM f"( x )( b ~ x ) 3 V f"(x)M(x) M(x)(b - x) 



dx 



v U(x) 

x£d(K b n[a,b]) 

U(x) 



+ A 2 (o) 



A', 



f"{x)M{xf 



i + v7"(x)m(x; 



1 + |M'(x) 



/"(x)M(x) 



f"(x) 2 M(x) 3 



1 + ^f"{x)M{x] 



and that A 2 (6) ; A 3 (&), and A^o) a// tend to /or some sequence ofb's tending to 
oo (The sequence may be different for diff eren t functions). The function A 2 (x) and 
the value of b are as defined in Theorem 4jJ_ In addition, assume that the sums 
and integrals in A3 (a), A 4 (x) ; and the K terms all converge as b tends to infinity in 



Theorem J±A 
Then 



9(n)e(f(n)) 



a<n<b 



■ E* 

f'(a)<r<f(b) 



g(xr)e(f(x r ) - rx r + §) 



-V(b) + c + A 



as a remains fixed and b tends to infinity, where c is some constant and 

A = 0(A 1 {b) + A 2 (6) + A' 3 (6) + A' 4 (6) + A 5 ) 
U(x) 



A f 



+ 



f"{x){x-af 
U(x) 



1 + 



+ 



f"(x)M(x 
1 + ^f"(x)M(x 



f"(x)(x-a) 

1 + \M'(x) 



f"(x)M(x) 



dx 



f"{x)M{xf 

+ K{J [b, 00), W , r ) + K{ J±[6, 00), W+, r+) + /C(J±[6, 00), W.,r_) 
g"(xf 



+ 



E 



g'{x)f"{xf 



and all other functions are as in Theorem 4-1 The sets Jo[b, 00), J±[b, 00), and 
J n uii[b, 00) all refer to the similarly named sets from Section 4-2 corresponding to the 
set [b, 00) instead of the set [a, b] . 

The implicit constants in the big-0 terms are dependent upon the constants G%, 
C2- , C4, Do, D\, D 2 , 6 in condition (M). 

An analogous result holds if b is fixed and a tends to —00. 



Much of the accuracy gleaned from Theorems 4.1 and 4.2 comes when f"(x) is 
small. If f"(x) is large, instead, then T>(x) shifts from being an explicit term to a 
big-0 term. However, with additional work it can be made somewhat explicit. 
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Theorem 4.4. Assume the conditions of Section 4-2 hold. Suppose that M (a), f"(a) > 
1 and let C and L be real numbers satisfying 

f"(a)- 1/2 < C < M(a) and Vf"( a ) < L < f"( a ) • mi 11 ! 1 , c }- 

t4/so let e = (a) and e' = (f'(a)). 
Then we may replace 



V°(a) = O U{a) 1 + — — + 



M (a) ^r{a)M(a) 

in the statement of Theorem 4J_ or \4-£\ with 
(17) 

'U{a)f"{a)C A L U(a)L U{a)f"{a) U{a) U{a) 



P » = P °< a > + °( M(a) + Mc + L 2 
where T> (a) can be estimated multiple different ways. 
• If e = 0, then 

'u(a)\e'\L , f/(a)| e '|(l + | e '|C)log(l + /») 



+ 



+ 



f"(a)C 2 M(a) 



Vn(a) = O 



+ 



+ U(a)f"(a)\e'\ 3 C 4 



f"{o) 
If e' = 0, then 

V (a) = i/j(a)g(a)e(f(a) - f(a)a). 
If e' = and \e\ > C, then 



V (a) = O 



U{a) 



{\e\-C)Jr{a) 



Suppose that e' = and M(a) < /"(a) 7 . We can then optimize the above results, 
so thatT>°(a) equals 

( ip(a)g(a)e(f(a) - f'(a)a) 

+ 0(U(a)f"(a)^\e\ 2 / 3 ) 



+ 



U(a) 



O 



U(a) 



U(a) 

M(a) 2 / 15 /"(a) 1/15 ' M{a) 



+ 



f"( a y/ 3 \e\ 2 / s 
+ 



U(a) 



+ 



U{a) 



< f"(a)-^ 



\e\ > f'(a)- 1 / 2 



M(a) 2 / 15 /"(a) 1/15 M(a) 

In all cases, the implicit constants depend on the constants associated to condition 
(M) and on the implicit constants in the bounds on C and L. This theorem holds 
true if a is everywhere replaced by b. 
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Remark 4.5. The results of Theorems 4J3 and \4-4\ can be combined to provide a 
more explicit T>(b) term as b tends to infinity. 

5. Necessary lemmas and useful propositions 

It is crucial to the van der Corput transform to have very good evaluations of 
integrals that contain stationary phase points. The most powerful results known to 
the author are those of Redouaby and Sargos in [30J . 

Lemma 5.1. (Lemma 2 in |30| ) 

Consider a function f : [a, b] — > R that is C 4 , and postive real numbers T, M, C 2 , 
C 2 -, C 3 , C 4 , Cm such that 

f"(z) > C~}TM- 2 fora<z<b 

f U) (z) < CjTM-i fora< z <b and j = 2, 3, 4 
b - a < C M M 

f'{c) = for some c in [a, b] 

Moreover, let 5 < 1 be positive, and define 

C2-C3 

IfO<c-a< (3M, then 

e(/(c) + l/8) / (3) (c)e(/(c)) e(/(o)) 



e(f(x)) dx 



1\[n^ 6vrzf' (c) 2 2m f (a) 

+ {f^w) +0 {f^) 



where 



1 fin \ , lM m ff(c + t)-f(c)\ 1/2 

a = -j (c) and ip(t) = sgn(t) ' 1 



2" w r w w V a 

IfO<b-c< f3M, then 

b dm) dx = <m±m + w±yw + <^ 



Z^fFic) Qmf"(c) 2 2mf'(a) 

r 4 \ (Ad 

T 2 (b-cf J +U \T 3 /\ 
The implicit constant is dependent upon the C 's and S. 

Remark 5.2. The restriction 5 < 1 in condition (M) guarantees that \if)(t)\ x \t\. 
See Lemma 4 in [30] . 

Also, the above result — which may be seen as a one-sided form of Huxley's result 
in [UJ — is a corollary to a more general result (Lemma 6 in [30] ), which makes the 
term 0{M 4 /T 2 \ji — c| 3 ) explicit. 
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Redouaby and Sargos include the additional assumption that T, M are both > 1 ; 
however, this condition is never used in the lemmas we cite, and so may be safely 
ignored^ 

Lemma 5.3. (Lemma 9 in |30j ) 

In addition to the hypothesis of the previous lemma, suppose that we have a function 
g : [a, b] — > C that is C 3 , and suppose there exist positive real numbers U, N, Do, D\, D2 
such that 

g ij) (z) < DjUN~ j fora<z <b and j = 0, 1, 2. 



Then 



<l[r)<UV))<h- <l(r) I ,(/(.,•)) 1,1 + ^Lf'la) G{f{a)) 

g'(c) , ^ ( UM 2 ( ^ t M 



'g(x)e(f(x)) dx = g(c) J e(f(x)) dx + ^~J^ </(>>)) 



This result still holds if we weaken the assumptions of the previous lemma so that 
f is only C 3 instead of C A , and the condition on f^\x) may likewise be ignored. 



Proposition 5.4. Suppose the hypotheses of Lemmas \5. 1\ and \5. 3\ hold and suppose 
a < c < b. Then 

g(x)e(f(x)) dx 

, , e(f(c) + 1/8) g(c)f^(c)e(f(c)) g{a)e{f{a)) 



and 



9{c) 
+ 



6mf"(c^ 2 



2mf'(a) 



2mf"(c 



7 /-e(f(c))+S(a), 



g(x)e(f(x)) dx 

, . e(/(c) + 1/8) , g(c)f^(c)e(f(c)) , g(b)e(f(b)) 



9{c) 



Qui f "(c) 2 



+ 



2m f '(b) 



12 



2m f "(c) 

This was confirmed in private correspondence with Sargos. 
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where 



S(y) = O 



UM 4 



T 2 \y 



+ 




Proof. We just apply Lemma 5. 3| followed by Lemma 5T 



□ 



The advantage of our condition (M) is now the following. 



Proposition 5.5. Suppose condition (M) holds, that f has a stationary phase point 
at x G J, and that x G [a, /3] C [x — M(x),x + M(x)\ R J. Then the conditions of 
Lemmas \5.1\ and 15.31 hold with 

C 3 = 1, M — -Mix), N = M(x), 

V 

C M = 2 V , U = U(x), T = r] 2 f"(x)M(x) 2 , 

and with CV, C2, C±, Dq, D\, D2 and 5 equal to the similarly named constants in 
condition (M). 

Morever, with these definitions, the error terms become 



£{a)+£{(3) = 



U(x) 



f'(x) 2 \x 



a>\ 



O 



U{x) 



f"(x) 2 \x- /3| 3 



O 



U{x) 



f(x) 3 / 2 M(xy 



where if x = a (resp., x = (3) then the first (resp., second) big-0 term disappears. 

However, we will also need good estimates on oscillatory integrals without a sta- 
tionary phase point. The first and second derivative tests below are well-known and 
often accompany any discussion of the van der Corput transform. Following them we 
give some more exact estimates that make use of the integration by parts technique 
we described in Section HI 

Lemma 5.6. (First derivative test — Lemma 5.1.2 in 

Let f{x) be real and differentiable on the open interval (a, (5) with f'(x) monotone 
and f'(x) > k > on (a, 0). Let g(x) be real, and let V be the total variation of g(x) 
on the closed interval [a, (3} plus the maximum modulus of g(x) on [a,/?]. Then 

,/3 

g(x)e(f(x)) dx < 



V 



7TK 



Lemma 5.7. (Second derivative test — Lemma 5.1.3 in [12] ) 

Let f(x) be real and twice differentiable on the open interval (a, (3) with f"(x) > 
A > on (a, (3). Let g(x) be real, and let V be the total variation of g(x) on the 
closed interval [a, (3} plus the maximum modulus of g(x) on [a,/3]. Then 



g(x)e(f(x)) dx 



< 



AV 

7^' 
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Proposition 5.8. (Condition (M) version of first and second derivative tests) 

Suppose condition (M) holds. If for some x G J, we have [a, (3} C I x := J D [x 
M(x),x + M{x)\, then 



g(y)e(f(y)) dx 



U(x) 



and 



g(y)e(f(y)) d y 



min xe[a! 0] \f'(x)\ 
U(x) 



< 



V 7 ? 



uniform implicit constant. 



Proof. We can bound V in Lemmas |5.6| and |5.7| by 

g(a) + g(/3) + f \g'{y)\ dy « U(x) + 2M(x) ■ ^ « U(x) 



M(x) 

using condition (M). Also by condition (M), we have that the minimum of f"(z) on 
I x is at worst ^> f"(x). □ 

In some cases, the first derivative test needs to be made explicit, for which we have 
the following result. 

Proposition 5.9. Let f G C 3 ([a,/3]), and g G C 2 ([a,/3]), and define h r (x) by 

_ (f'( x )- r )g'(x)-g(x)f"(x) 
r[X) (f'(x) - r)* 

Suppose that f\x) ^ r on an interval [a, j3], and let 

K r (a,P) :=^2\hr(x)\ 

where the sum ranges over all x G [a,/3] where h' r (x) = 0. Then we have 



g(x)e(f(x) — rx) dx 



1/3 



-e(f(x) — rx) 



2iri(f'(x) -r) 
+ 0(K r (a,f3)) + 0(\h r (a)\) + 0(\h r (/3)\) 

Proof. We begin by applying integration by parts twice to our original integral. 

g(x)e(f(x) — rx) dx 



(18) 



g(x) 



-e(f(x) — rx) 



2m(f'(x) - r) 
V 1 (f(x)-r)g'(x)-g(x)f"(x) 



2m 



U'{x)-tY 



e(f(x) — rx) dx 
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(19) 



(20) 



9{x) 
2m{f'{x) 

1 



-e(f(x) — rx) 



_(2tt«) 2 



+ 



h r (x)e(f(x 
d 



rx) 



dx 



h r (x) 



dx 



The two terms on line (19) are bounded trivially by 0(h r (a)) + 0{h r {b)). 

The remaining integral on line (20) is bounded by the total variation of h r (x), 
which is at most the modulus of h r (x) at every critical point added to |/i r (a)| and 
\h r (/3)\. K r (a,{3) is, by definition, the sum of the moduli h r (x) at every critical point 
(since the derivative — by the assumption that f'(x) ^ron [a, (3} — always exists), so 
the proof is complete. □ 

By making use of the following proposition, we can sometimes avoid the K r term 
in Proposition 5.9 An alternate form of the resulting error terms can be found in 
Lemma 5.5.5 of 



Proposition 5.10. Suppose condition (M) holds. If e— ±1, then for all x such that 
x, x + e • M(x) G J and all r such that 

f'(x) — r = or sgn(f'(x) — r) = sgn(k), 

we have 

\f'{x + e • M{x)) - r\ > M{x)f"(x). 
Proof. By Taylor's remainder theorem, we can write 

f\x + e • M{x)) -r = (f(x) - r) + e • M{x)f"{x) 

1 



(21) 

with implicit constant 1. 
By condition (M), we have 



+ -M{x) 



max f^ 3 \y) 

yE[x,x+kM(x)] 



max 

yd[x ,x+e- M (x)] 



\f i3) (y)\ < 



r]f"{x) 
M(x) ' 



Therefore the big-0 term in (21) is at most r)M(x)f"(x)/2. Since rj < 2 and since 
f'(x) — r is either or shares the same sign as e, this gives the result. □ 

Proposition 5.11. Suppose f and g satisfy condition (M) for some function M(x) 
and its associated constants. Let e = ±1. Assume that x r is not in the interval from 
x to x + e • M(x), and that x + e • M(x) is farther from x r than x is. Then we have 
that 

g{y) 



'x+e-M(x) 



g(y)e(f(y) - n/) dy 



2iri{f'{y) - r 



■e(f(y) - r y) 



y=x+e-M(x) 
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U(x)f"(x) 



+ 



|/'(z)-7f 



o 



U(x) 



We also have the following bound: 

g(y) 



2ni(f'(y) - r 



e(f(y) -ry) 



O 



M{x){f'{x) -r) 2 
U(x) 



M(x)f"(x) 



y=x+e-M(x) 

Proof. As in Proposition |5.9[ we apply integration by parts twice to our starting 
integral 



g(y)e(f(y) - ry) dy 



x+e-M(x) 



9(y) 



~Af{y) - r v) 



y=x+e-M(x) 



(22) 

The relation 



2m(f'(y) - r 
+ 0(h r (x)) + 0(h r (x + e • M{x))) 

d (f(y)-r)g , (y)-g(y)f"(y) 



O 



x+e-M(x) 



dy 



{2-KiY{f'{y)-rf 



dy 



g{y) 



2m(f'(y)-r 
holds by Proposition |5.10[ 



-e(f(y) - r y) 



O 



y=x+e-M(x) 



g{x) 



M(x)f"(x) 



Using condition (M), we have the following additional bounds. 



h r (x) = O 
h r (x + t ■ M(x)) = O 



g( x ) 



o 



g{x)f"{x) 

\f'{x) — r| 3 



M(x)(f (x) -r) 2 

g'(x + e- M(x)) 
(f'(x + e-M(x)) -r) 2 

g(x + e ■ M(x))f"(x + e ■ M{x)) 



O 



(f'(x + e- M(x)) -r) 3 



O 



g(x) 



o 



g(x)f"(x) 

\f'(x)-rf 



M(x)(f'(x) -r)\ 

The last equality holds since \f'(x + e- M(x)) — r\ is bigger than \f'(x) — r\ by our 
assumption about the relative placement of x, x r , and x + e ■ M(x): for example, if 
x r < x < x + M(x), then r < f'(x) < f'(x + M(x)). 



Finally, to estimate the integral in (22) first notice that 
d (f'(y)-r)g'(y)-g(y)f"(y) 



dy 



(2my(f'(y) - r) 3 



O 



g"(y) 
(f'(y) - r) 2 



o 



( g'{y)f'{y) 
W'(y)-r) 3 



, n f g(y)f i3) (y) \ , n ( g(y)f"(y) 2 
W'iy)-r) a J W'(y)-rr 
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On the interval between x and x + e • M(x), the maximum of g" (y) / (f (y) — r) 2 is 

g( x ) 

< M(x) 2 (f'(x) -r) 2 ' 

and since we are integrating over an interval of length M(x), the contribution of this 
term is at most 



M(x)(f'(x)-r) 



Similarly, one can show that 

g'(y)f"(y) 



o 



x+e-M(x) \\f'(y)-r\ 3 



+ 



g(y)f (3) (y) \ 
\f'(y)-r\ 3 J 



dy = 



g(x)f"(x) 



\f>(x)-r\* 



We estimate the remaining integral by using f"(y) = 0(f"(x)): 



I 



g{y)f"{y) 

x+e-M(x) (f'(y) ~ r) 



i dy = 0[ g(x)f"(x) 



O { g(x)f"(x) 

g{x)f"{x / 



f"(y) 



x+e-M(x) (f'(y) - rY 
1 

(/'(•'/) - , ') :i J J 



dy 



O 



\f'(x) — r| 3 / 

Remark 5.12. If g(x) is a constant, then we may remove the term 

U(x) 



□ 



O 



M(x)(f'(x)-r) 



from the statement of Proposition 5.11, This would remove the term U(x)/M(x)\\f'(x)\ 
from A2(x) in Theorem 4-1 ■ However, very often one finds a term of this same size 
occuring in A%(x) regardless. 

The primary technique we shall use to evaluate sums will be a variant the Euler- 
Maclaurin summation formula; however, the Euler-Maclaurin formula itself sums over 
all integers in an interval [a, b], while we will often want to sum over all values x in 
some interval [a, (3] for which F(x) is an integer (for some function F). 

Lemma 5.13. (Euler-Maclaurin summation, first derivative version — page 10 in 
Suppose f is a differentiable function on [a, (3] then 

= / dx+ s(x)f'(x) dx 



a<n<l3 
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Proposition 5.14. Suppose G is a real-valued, differentiable function on [a,b]. Sup- 
pose F is a differentiable function such that F([a,b}) = [a,/3], then 

dy 



G(F-\n))= I G(y).\F'(y)\dy + o( [ \G'(y)\ 
+ O (\s(a) ■ G(a)\) + O (\s(b) ■ G(b)\) 



a<n</3 



+ Yl H*)-G(x)\ 

\x is a local extremum of F on [a,b] 

where G , (F _1 (n)) is a sum over all G(x) for x G 

Proof. We apply the first derivative version of Euler-Maclaurin summation to each 
interval where F is monotonic (and hence, F _1 is 1-1). Then we apply a change of 
variables with y = F^ 1 {x). The sum 

of E N*)-o(*)l 

\ x is a local extremum of F on [a,b] 

bounds the contributions of s(a)f(a) and —s(/3)f(/3) arising from the endpoints of 
the intervals of monotonicity. □ 

Lemma 5.15. (Lemma 5.4-2 in [12p*| — Poisson summation) 
Let f G C 2 ([a, b}) for real numbers a < b, then 

R fb 

f{ n ) — hm / f(x)e(rx) dx. 

a<n<b C ° r =-R^ a 

Proof. We rewrite the Euler-Maclaurin summation formula as 

f(n)= f f{x)dx-iP{/3)f(/3)+iP(a)f{a) + f ^{x)f{x) dx. 

a<n<p Ja Jol 

But we also have 

" ' — sin(2irrx) 

/ (x) dx 



rP ] rP " 

/ ib(x)f'(x) dx = lim / > 

Ja 7T -R— s>oo J a ^ 

R „/3 

= ij)(/3)f(a) — ijj(a)f(a) + lim N / 2cos(2nrx)f(x) dx 

R .p 

— — ^( a )/( a ) + li m / / (e(rx) + e(—rx))f(x) dx, 

J a 

13 Huxley assumes a and b are integers, so we provide the proof when a, b are real. 
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by integration by parts, which completes the proof. □ 

We end this section with two more specific propositions. The first will give us good 
bounds on the partial sums and tails of ip(x,e), our modified sawtooth function, as 
well as guarantee its convergence. The second will show that the integrals in Theorem 



4.1 involving the W and r functions will converge. 



Proposition 5.16. Suppose (3 > a > 1, e E [—1/2, 1/2], and j;6l Then 

\- e(mx) ( . J 1 

> = O mm < 1, — 77— r . 

m + e \ Oi\\x\\ 

a<\m+e\</3 V 

where the implicit constant is uniform. Also, the sum is convergent as (3 tends to oo 
and the same bounds hold in this case. 

Proof. Without loss of generality, we may replace that x with (x) and so may assume 
that x E [—1/2, 1/2). We start with the assumption that (3 is finite. 

Next, we remove all appearances of the e from the sum. We use C a , x ^ and C' a x e 
to denote constants which only depend on the variables a, x, and e. 



E ^= E ^ + c w + ofi 

L — ' m + e L — ' m + e \3 

a<\m+e\<f3 o<M</3 



E 



e(mx) \ ^ 

eimx) 



m v ' V m(m + e) 

a<|m|</3 o<|m|</3 



E ^-E^fV^Vofi 

a<\m\<p a<\m\ \ \ < J / \r 



+ C a .x.e + O 



a<|m|< 

Here, Cq, jXi(E contains the terms where \m\ < a < \m + e| or \m + e\ < a < \m\. In 
particular, it is no larger than 2/ (a — 1/2) = 0(1/ 'a). Similarly, C' axe is dominated 
by the sum of terms 1/m 2 , and hence is also 0(1 /a). 

Now we pause a moment to show that we may let (3 = oo with no problems of 
convergence. The sum 

% ^ e(mx) 
' m 

l<|m| 
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is, up to a constant multiplier, the Fourier series for the sawtooth function ip(x) and 
converges for all x. This implies that 



hm y 



e\mx] 



>oo — ' m 

a<\m\</3 



converges for all a and x. Therefore, since 

I ^ ejmx) v e(mx) , , 

\a<\m+e\<l3 ot<\m\</3 



we have that 



converges for all a, x, and e. 
If x = then 



and hence 



e{mx) 
- / -— ' m + e 

a<|m| 



y- e ( ma; ) = 

a<M</3 



y- e H = o 
^-^ m + e \a 

a<\m+e\</3 v 

So we may assume for the remainder of the proof that x 7^ 0. 

Now we remove the absolute value in the condition on the sum and apply summa- 
tion by parts. 

y-^ e{mx) 2. y^ sin(27rmx) 
^ m - / -— ' m 

a<|m|</3 a<m</3 



2 > E Marfa)) 

a<m</3 v 7 \a<fc<m / 

El sin(7r(m + [a] )rr) sin(7r(m — [a] )x) 



m 2 + m sin(7nr) 

Since we assumed x G [—1/2, 1/2], we have that sin(7r:r) x x. For the numerator, we 
use that |sin(y)| < min{l, \y\}. This gives the following estimate, which completes 
the proof. 



y- e ( mx ) = ( y- J_ 



e(mx) I 1 . f 1 2| . . 

2 mm ^ 71 — m ' m \\ x \ 



m \ <■ — ' I \\x 

a<\m\</3 \a<m</3 ^ 11 
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o E w +o E 

\ a<m<l/||x|| / \ m>max{a,l/||a;||} 

'0(1) ifa||z||<l 

^ I ~~ n — rr I otherwise 

ax 



m 2 ||x| 



□ 



Proposition 5.17. Suppose the conditions of Section \4-£\ hold, then the function 
l^o( x )l + l^o^) ' r o(^)| is integrable on J and the function \W+(x)\ + |Wi(a;)| + 
|W+(x) ■ r',(x)\ + |W_(x) • r'_(x)\ is integrable on J±. 



Proof. Recall that 



H = gf® + 3g'f" 
G = Ugg"{f"f 

{2g"fg' (2g") 3 (f"g) 



W± 



r± = f 
ro = f" 



(H ± VH 2 - G) 2 (H ± VH 2 - GY 
H 2 f^ 

27g(f") 5 

H ± VH 2 - G 



2g" 

un 2 



H ' 

The functions H and G are continuous and bounded on J. The functions Wo and r 
are continuous inside Jo, and the functions W± and r± are continuous inside J±. The 
only possible barrier to integrability is the presence of a zero in the denominator of 
the function at an endpoint of an interval. 

The only terms in the denominator of Wq are g and /", which are both bounded 
away from on J . The only terms in the denominator of Wq are g and /" again, 
but r' has a factor of H 2 in the denominator; however, in |Wo(a;) • r (x)\, the H 2 in 
the denominator of r' Q is cancelled by the factor of H 2 in the numerator of Wq(x), so 
this too has no zeroes in the denominator at the endpoints of Jo- 

The case when we are considering an interval (a', b') where H 2 (x) — G(x) = but 
H(x) does not equal nor tend to is analogous. So we will assume we are in the case 
where H 2 (x) — G(x) does not equal nor tend to on the interval in question for the 
remainder of this proof. 

We have 



r'±{x) = f" 



H ± VH 2 - G 
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and 



H± VH 2 -G\~ 4 fH±VH 2 -G\' /c „ . 

(/ 9) 



In addition 



2^" J V 2(7" 
} H±VJP~G y i ( H±y/H*-G \ , 
" A 2^" J v 2g" ]IJ - 

H±y/H 2 -G\' H±VH 2 -G 



(23) I I = -g® 

1 W J 2g" 2 



2g"VH 2 - G \ \ J 2 J 



So, the only terms in the denominator of r' ± are g" and \JH 2 — G, and the only terms 
in the denominator of W' ± are H ± \/H 2 — G and \/H 2 — G as all the g" terms cancel. 

The technique required will change slightly depending on whether we consider the 
+ terms or the — terms. As G tends to (which is equivalent to g" tendi ng t o 0), 



H ± V H 2 — G will tend to either 2H or 0. By the assumptions of Section 4.2 , if G 
tends to at the endpoint of an interval in J±, then H cannot tend to at the same 
point. 

If the sign is chosen so H ± V H 2 — G tends to 2H, then the terms H ± \/H 2 — G 
and V ' H 2 — G in the denominator of W'± (x) do not vanish, and the two copies of g" 
in the denominator of r' ± (x) are cancelled by the two copies of g" in the numerator 
of W±(x). Therefore, in this case, there are no zeroes in the denominator. 

For the remaining case, when H ± \/H 2 — G — > 0, we have that 



H ± VH 2 - G 




2g>- 

(24) 



H 

where the implcit constant depends on the size of H near this point. In particular, 



since g and /" do not approach at the endpoints of J±, equation (24) implies that a 
copy of H ± y/H 2 — G tending to in the denominator of a function can be cancelled 
by a copy of g" in the numerator to prevent the presence of a zero in the denominator. 
In particular, W±(x) will not have a zero in the denominator. 
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By (24), the only term in |W4(x)| + |W±(x) • r' ± (x)\ that could produce a zero in 
the denominator is the derivative 

H± VH 2 -G\' 



(25) 



2<7'' 



Therefore, it suffices to show that (25) has no zeroes in the denominator at an end- 
point of J±. 

By applying (24) to line (23) and noting that in this case ±(H 2 — G)~ 1 / 2 = 
-H- l (l + 0(G)), we have 

H ± VH 2 -G\' 3g(f") 2 gW ( \ 



2g» 



9" 



± 



3g(f") 2 H' 

H y/H 2 -G) ' 2VH 2 - G 



± 12g'(f") 2 + 2Agf"f^ + 



AVH 2 - G 



9" 



■O 



=F 



2VH 2 - G 
± l2g'{f") 2 + 2Agf"f^ 

Ay/H 2 - G 
3g(f") 2 H> l2gi{f"Y + 2Agf»fW 



2VH 2 - G 



± 



A^H 2 - G 



+ 0(1). 



Since \/H 2 — G does not tend to on 4, the derivative (25) has no zeroes in the 
denominator. □ 

6. Proof of Theorem 14.11 



6.1. The initial step. We begin by applying Poisson summation (Lemma 5.15). 
This gives 

* R r b 

} \ 9{n)e{f{n)) = lim V / g(x)e(f(x) - rx) dx. 

a<n<b r =-R Ja 

We wish to alter the end-points of some of these integrals. Define a r , the new left 
endpoint, by 

( a + M(a) if < /'(a) - r < f"(a) 
(26) a r := la- M(a) if > /'(a) - r > -f"{a) 

I a otherwise 

and b r , the new right endpoint, by 

'b + M{b) if < f{b) - r < f"(b) 
b - M{b) if > /'(&) - r > -f"(b) . 
b otherwise 
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Note that f'{a) — r is positive, if and only if x r lies to the left of a and vice- versa. 
The transformation from [a, b] to [a r , b r ] has the effect that if x r is close to — but not 
equal to — an endpoint a or b, then we shift that endpoint away from x r by a distance 
M(a) or M(b), respectively. 

In the statement of Theorem |4.1 m x is defined by 

m x := |Z n {f{x) - f'(x), f\x) + f\x)) \ {f(x)}\ . 

The values m a and m& thus count the number of r such that a r 7^ a and b r 7^ b, 
respectively. 

After altering the end-points, we have 

R rbr 



9(n)e(f(n)) = lim / g(x)e(f(x) - rx) dx 

D J 0, r 



a<n<b r=—R 



(27) — Jim I y + ^ J g{x)e{f{x) — rx) dx 



-R 



br 



Since there are only finitely many values of r for which a r ^ a and 6 r 7^ 6, the sum 
on line (|27|) is finite and so we may omit the limit and sum from —00 to 00. 



First, consider those integrals in the sum on line (27) arising from < \f'(a)—r\ < 
1. If (f'(a)) = or if m a = 0, there are no such integrals, and so the total contribution 
of this (null) set of integrals is 0. Otherwise, there are at most 2 such integrals, each 
of which by Proposition 5.8 (our variant of the first and second derivative tests) is 
bounded by 



(28) O min 



U(a) U(a) 



vTW ll/'WII 



Thus the contribution of these integrals is 0(Ai(a)). 

There can only be additional integrals in line (27) if f"(a) > 1, so we will assume 
as such when bounding them. 

Consider next those integrals from a r to a arising from 1 < \f'(a) — r\ < a/ f"{p)- 




There are a/ f"(a) such integrals, each of which, by Proposition 5.8 is bounded by 

o( u{a > 

So the total contribution of these integrals is bounded by 0(g(a)). 

Lastly, we consider the sum of those integrals in line (27) arising from f"(a) < 
\f(a) — r\ < f"(a). In particular, let 

pa 

Si := ^2 / 9(x)e(f(x) - rx) dx 



7"W<l/'W-r|</"(") 
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denote the sum in question. We apply Proposition |5.11| to each integral in Si to 
obtain 



(29) 
(30) 
(31) 



9W 



- 27ri(f'(a) - r) 
U(a 



e(f(a) - ra) 



r 



f"(a)M(a 
U(a 



M(a)(/'( a )-r) 2 



O 



U(a)f"(a) 
(/'(a)-r)» 



Here, each sum is over all r satisfying ■sj f"(a) < \f'(a) — r\ < f"(a)- 

For the sum on line (29), we apply a change of variables m + e = /'( a ) — r > 
where e = (/'(a)), and then apply Proposition 5.16 The sum is then bounded by 
0(U (a) I \J f" (a)) = 0(U(a)), since we are under the temporary assumption that 
/»>!■ 

The sum on line (30) has at most 2 • f"(a) terms, so is bounded by 

U(a) 



O 



M(a) 



By Euler-Maclaurin summation (Lemma 5.13), the sum on line (31) is bounded by 

O I U(a) 



1 



Jna)M{a) 



+ 1 



Thus the total contribution of the integrals from a r to a in line (27) arising from 
1 < |/'( a ) - r l < f"( a ) is bounded by 



O U{a) 1 + 



1 



+ 



1 



where 



M(a) yTf^Mia) 
0(%i(«)), 



^>i(a) = U(a)[l + — + 



Similarly the contribution of the integrals from b r to b in line (|27j) is bounded by 

0(V° fl ^(b)) + 0(Ai(b)). 

Therefore, 

R rb r 

y~) g( n ) e (f( n )) = i im y~] / 5'(a ; ) e (/(a ; ) - tx) dx 



a<n<b 



(32) 



r=-R " a r 

+ 0(P>> 1 (a))+0(^„> 1 (6)) 
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+ 0(A 1 (a) + A 1 (b)). 

6.2. The half stationary phase estimate. Suppose that f'(a) = r£Z. We will 
need to estimate the integral at this particular r in (32) separately. 

We shall break the integral into a piece surrounding the stationary phase point and 
an integral away from the stationary phase point. In particular, we write 



g(x)e(f(x) — rx) dx 



+ 



g(x)e(f(x) — rx) dx, 



where (3 r = min{a + M(a), b r }. Since x r = a is to the left of b, we have that b r either 
equals b or b + M{b). 

By Propositions 5.4 and 5.5, we have 

g(y)e(f(y) - ry) dy 

g(a)e(f(a) -ra + 1/8) , g(f3 r )e(f((3 r ) - r(3 r ) 



+ 



2^r[a) 27u(f(/3 r )-r 
g(a)f^(a)e(f(a)) g'(a)e(f(a)) 



+ 



Qmf"(a) 2 
U(a) 



+ 



2%if"(a) 
U(a) 



+ 



U{a) 



f"(a) 2 (b-af 



(33) 



J"{a) 2 M{af f"{afl 2 M{a) 2 , 
g(a)e(f(a) -ra + 1/8) , g(0 r )e(f(Pr) ~ r(3 r 



+ 0(A 1 (a)+A 2 (a)). 
Here we used the bound 



+ 



2m(f'(/3 r ) - r) 



+ V*(a) 



1 



(A- - a) 3 



O 



1 



(b-a) 



O 



Af(a) 



Similarly, if /'(&) = r G Z, then we dissect the integral from a r to 6 in (32) in a 
similar way, obtaining, 



+ 



g(x)e(f(x) — rx) dx 



(34) 



g(b)e(f(b) - rb + 1/8) ^(a r )e(/(a r ) - ra r ) 



-£>*(&) 



2/p(6) 27Tz(/' («r) - r) 

+ 0(Ai(6)) + 0(A 2 (6)) + / ' g(x)e(f(x) - rx) dx, 

J a r 

where a r = max{6 — M(6), a r }. 
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6.3. The full stationary phase estimates. For all remaining integrals with a 
stationary phase point in (32), we may assume that \f'(a) — r\, \f'(b) — r\ > 0. As 
before, we denote the stationary phase point corresponding to a given r by x r . For 
each such r 6 (f'(a), f'(b)), we write the corresponding integral in line (32) as 

"b r 



g(x)e(f(x) — rx) dx 



+ 



+ 



Pr 



g(x)e(f(x) — rx) dx, 



where a r = max{x r — M(x r ),a r } and /3 r = min{x r + M(x r ),b r }. By construction, 
[a r ,/3 r ] = I Xr PI [a r ,6 r ] with I x as in the statement of condition (M). 
For this section, we shall focus on the contribution of the middle terms. Consider 



f'(a)<r<f'(b) 



Pr 



g(x)e(f(x) — rx) dx. 



Applying Propositions |5.4| and |5.5| with c = x r to S 2 , this becomes 
(35) S 2 = 



/'(a)<r</'(6) 



g(x r )e(f(x r ) - rx r + |) 



(36) 
(37) 
(38) 



/'(a)<r</'(6) 



gr(X r 



^ V/"(a; r ) 2 V(^-«r) 3 



+ 



f(a)<r<f'(b) 



g{x r 



- E 



/'(a)<r</'(6) 



f'(x r ) 2 M(x r ) 3 

g(a r )e(f(a r ) - ra r ) 
2m(f'(a r ) - r) 



1 + y/f"(x r )M(x r )^ 

9(Pr)e(fWr) - r(3 r ) 



+ 



E 



/'(a)<r</'(6) 



2m(f'((3 r ) - r) 



The sum in line (35) added to the first terms in lines (33) and (34) — if they exist— 
sum to 

9(x-r)e(f(x r ) - rx r + |) 

f'(a)<r<f'(b) Vf'( X r) 

the main term in the van der Corput transform. 

To begin evaluating the sum in line (36), consider values r such that < x r — a < 
minjC^ 1 , M(a)}. By the mean value theorem, we can write x r — a = (r— /'(a))//"(C) 
for some ( G (a,x r ) C (a, a + M(a)). By the bounds of condition (M), the r under 
consideration must also satisfy < r — f'(a) < f"{a), which, in turn implies that 
x r — a r = x r — a + M(a) > M(a). 

Therefore, we split the sum over 

g(x r ) 



f {x r ^{x r Oj r )^ 



48 



J. VANDEHEY 



in line (36 ) into two pieces. The first piece, over all r such that x r — a < min{C 2 1 , M(a)}, 
has at most m a terms, with each term bounded by 

m ?M uio) 

K ' f"{x r ) 2 M{af f"(a) 2 M(a) r 

Thus the total contribution of these terms is at most 

U(a)m a 



O 



0(A 2 (a)), 



J"(a) 2 M(a)\ 

since m(a) = 0(1 + /" (a)). 

The second piece, over all remaining r between f'(a) and f'(b), is bounded by 



(40) 



fii 



a<x r <b 



f"{x r )\x r -af 



where a is the smallest value in the interval [a + min{M(a), C^ 1 }, b) such that f'(a) 
is an integer. Note that if no such integer exists, then this sum is empty. We now 
apply Proposition 5.14 to (40) with 

G{y) 



g(y) 



and F(y) = f'{y) on the interval [/'(a), fib)]. Using the bounds from condition (M) 



liberally, we see that (40) is bounded by 

r-b 



(41) 
(42) 



O 



U(x) 



f"(x) 2 (x-a) 3 

U(a 



f"(x) + 



O 



f"{a) 2 {a-af 



M(x) x — a 

+0 , m 



dx 



f"(b) 2 (b-af 



0(A 3 (a)). 



Since f"(x) > on [a, b], f has no local extrema on this interval, so the final term 
from Proposition 5. 14| does not appear. 

We have a similar bound for the second sum in line (36): 

g(x r ) 



E 



(f"(x r )) 2 (b r - x r f 

f'(a)<r<f'(b) W V T " V rJ 



0(A 2 (6)) + 0(A 3 (6)). 



(43) 



(44) 



We apply Proposition |5.14| again to the sum on line ( |37| ) to get 

1 + Vf"(Xr)M(x r ) 



g{x r ) 



o 



1 + 



1 + \M'(x)\ 
f"(x)M(x) 



dx 
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+ 



< 45 > +0 (?^w( 1 + ^^)) 



J"(b) 2 M{bf 
= 0(A 4 ) + 0(A 2 (a) + A 2 (6)). 

Now consider 
S3 : = lim y / g(x)e(f(x) — rx) dx 

\r\<R Jar 
r</'(a) or r>/'(6) 

/'(a)<r</'(6) v w VK ' 7 17 Pr 

+ e or^ww-")*- ^:;:^ )- 

f'(a)<r<f'(b) yJar \J \ rj / / 

Note that the remaining terms of S2 in line (38) as well as the second terms from 
lines (33) and (34) — if they exist — appear in S3. 
We have thus far shown that 

v^* ( , (f( u v-* g(x r )e(f(x r ) -rx r + D t»*^^t»v ^ 

2^ g{n)e{f{n)) = ^ /7 — ^ — + £3 - 2? (6) + V (a) 



a<n<b f'(a)<r<f'(b) 

4 

+ O^^a)) + 0(V° r ^(b)) + ^ 0(Ai(a) + A<(&)) + 0(A 4 



j=i 



6.4. The remaining integrals. For each integral in 5*3, we apply Proposition |5.9 
This gives 

(47) + lim V O (h r (a r )) + lim V O (h r (b r )) 

r£f'(a) r^/'(6) 

(48) + J] (0(/ ir (a r .)) + 0(/ ir (/3 r ))) 

/»<r</'(6) 

+ Yl (0(K r (a r ,a r )) + 0(K r (f3 r ,b r ))) 

/'(a)<r</'(6) 

+ lim > 0(K r (a r ,b r )) 

R-¥oo — ' 

r<f'(a) or r>f'(b) 

where in each sum it is assumed that Irl < R. 
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We start estimating the first sum of (46) and write 

g(a r )e(f(a r ) - ra r ) 



— lim > 



2m(f'(a r ) - r) 



(49) 
(50) 
(51) 



E 

0<r-/'(o)</"(o) 



+ E o 

0<-(r-/'(a))</"(a) 

a(a)e(/(a)) 



p(a-M(a)) 
/'(a - M(a)) - r, 

g(a + M(a)) 



f'(a + M(a))-r 



2m 



E 



llm / en 

|r|<R ^ V 7 
|/'(a)-r[>/"(a) 



-ra 



There are at most m a terms on lines (|49j) and ( 50 ) , each with size 

o( m 

\f"(a)M(a) 
by Proposition 5.10 for a total contribution of 

U(a)m a 



O 



/"(a)M(a) 



C(A 2 (a)). 



We can evaluate the sum on line (51) explicitly in some cases. Suppose /"(a) < 
1 - ||/' (a)||. If f"{a) < ||/'(a)||, then we pull out the term where f'(a)-r= (/'(a)}, 
so that the remaining terms on line (51) form a modified sawtooth function. In 
particular, the sum on line (51) becomes 



f a(a)e(/(a) + [/'(a)]a)x 
1 



x 



+ V(a, (/'(a))) 



/»< ||/'(«)|| 



2m(f'(a)} 

j(a)e(f(a) + [/'(a)] a)^(a, (/'(a))) ||/'(o)|| < /"(a) < 1 - ||/'(a)ll 

which equals V°(a), when f"(a) < 1 - ||/'(a)||. Otherwise, if f"(a) > 1 - ||/'(a)||, we 
may bound the sum on line (51) using Proposition 5.16 by 0(U(a)/ f"(a)) = 0{U(a)). 
Recall that T>°(a) is not explicit when f"(a) > 1 — ||/'( a ) II; an d instead is just a big-0 
approximation with one of its terms of the form 0(U(a)). 

By a similar argument, one can show that all the sums on line (46) are equal to 

-V°(b) + X?» + 0(A 2 (a) + A 2 (6)). 



We break the sum from (47) into three smaller sums like so 
(52) lim V 0(h r (a r )) = lim V OfhJa)) 



[r|<fl 



k|<-R 
l/'(a)-r|>/"(a) 
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(53) 
(54) 



+ 0(h r (a - M(o))) 

0<r-f'(a)<f"(a) 

+ 0(h r (a + M(a)). 



0<f'(a)-r<f"{a) 



The sum on the right-hand side of line (52) is 



lim 

R— >oo 



E 



|r[<R 
\f'(a)-r\>f"(a) 



(f(x)-r)g'(x)-g(x)f"(x) 



(55) 



<C lim 



E 



|r|<fl 
l/'(a)-r|>/"(a) 



(f'(x) - r)3 
U(a) 



M(a)(f'(a)-ry 



U(a)f"(a) 



(/'(a)-r)» 



If ||/'( a )|| — or m a > 1, then the smallest \f'(a) — r\ could be in this sum is 
min{l/2, f"(a)}; otherwise, we will have a term where \ f'{a) — r\ = ||/'(a)||. Thus 
this sum is bounded by 



f U(a) 



M(a) 



min < 1 



1 



/"(«) 



< < 



U(a) 



+ U(a) min<j f"(a), — — 



U(a)f"(a) 



[M(a)\\f'(aW \\f'( 



if ||/'(o)|| = 0or m a >l 



otherwise 



< A 2 (a). 



By using condition (M) and Proposition 5.10, each term on lines (53) and (54) has 
size at most 



U(a) 



/"(a) 2 M(a) 3 ' 

and there are at most m a = 0(1 + f'(a)) such terms for a total contribution of at 



most 



O 



U{a) 



f"{a) 2 M(af 



(! + /») =0(A 2 (a)). 



By a similar argument both the sums on line (47) are bounded by 0(A2(a)+A2(6)). 

For the first sum on line (48), we have a r equals a r or x r — M(x r ). The sum of 
terms with a r = a r are bounded by the sum in line (47), which we just showed has 
size 0(A 2 (a)), so we only need to estimate the terms where a r = x r — M(x r ). By 
Proposition 5.10 and the bounds of condition (M), each such term is bounded by 



g(x r 



f»(x r yM(x r y 
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We may therefore bound the sum over those r with a r = x r — M(x r ) by a sum of 
terms g(x r )/f"(x r ) 2 M(x r ) 3 for all r e [/'(a), /'(&)] and obtain 

= (~)( A ) 



by the same argument as in line (43). 
Therefore, we have shown that 

\* \* q(x r )e( f (x r ) — TX r + I) 

9{n)e{f{n)) = £ ^ W /Zr^ ~ " P ^ + 

a<n<6 f'(a)<r<f'(b) V * ^- X ' r ' 

3 

+ J]0(A l (a) + A,(6)) + 0(A 4 ) 
i=i 

+ £ (0(# r (a r , a r )) + 0(# r (A, 6 r ))) 

f'(a)<r<f'(b) 

+ 0{K r (a r ,b r )). 

r<f'(a) or r>/'(6) 

The proof will therefore be finished when we show that the sum of all the K r terms 
is bounded by 0(A 4 ). 



6.5. Bounding the variation. The K r 's are the sum of terms of the size 

(/'(*) - r)g'{x) - g(x)f(x) 



hJx) 



(f'(x) - r)3 



at points x where the derivative with respect to x vanishes. We may safely ignore 
points where the derivative does not exist, since all the intervals that give rise to the 
K r terms do not contain stationary phase points. 

The derivative of h r (x) is (using f' r as shorthand for f'(x) — r) 

h , (x) = -gfrf (3) +9" (fr) 2 + Mf'T - yfrf" a" ifrf ~ # f'r + MfT 



"r V 



We set the numerator equal to and solve for f' r . 

First, suppose g(x) = 0. If g"(x) also equals 0, then the numerator is if and only 
if g'(x) also equals (f"(x) is never 0, and we assumed f'(x) ^ at all points in 
consideration); but if g(x) and g'(x) equal 0, then h r (x) also equals 0, so these points 
contribute nothing to K r . Thus the only contribution from points x where g(x) = 
come when g"(x) ^ and g'{x) ^ 0, i.e., those points in J nu n as defined in Section 



4.2 At these points, solving for f' r yields 3g'(x)f"(x)/g"(x), so the contribution to 
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the K r terms from these points is bounded by. 



« T 9 [x) 



For the remainder of this section, we will suppose then that g(x) ^ 0. There is no 
way to make the numerator equal if g"(x) = H{x) = 0, since both g(x) and f"(x) 
are assumed to be non-zero, or if H(x) 2 — G(x) < 0, since all the above functions are 
real-valued; otherwise, the solution is given by 



H ± VH 2 - G j 

(56) /; = < ur/ X " ! • 



where J± and Jo are as defined in Section 4.2 Plugging this value of f' r into h r , we 
obtain 



h r (x) 



W±(x) x E J± 

W Q (x) x e J 



where 



w ± = < 2 f )V {29 " )Hf " 9) - d 

(H ± VH 2 - G) 2 (H ± y/H 2 — G) 3 



W 



H 2 g> H 3 f"g 



(3g(f") 2 ) 2 (Uf") 2 T 

#2/(3) 



27s(f") 5 ' 

If we ignore the constraint that r be an integer for the moment, then we can imagine 
that (56) determines functions 



H±s/H 2 ~G 
r±{x) = f — on J ± and 

ro(x) = f — — on J . 

The contribution of the K r terms is then at most 

r~ 1 (n)£.J± rZ 1 (n)GJ± 
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We now apply Proposition 5.14 to see that the contribution of the K r terms is 
bounded by 

« K{J , W , r„) + £(J±, W+, r+) + K{ J±, W.,r_) 
<0(A 4 ). 



This completes the proof of Theorem 4.1 



7. Proof of Theorem 14.21 



We follow the proof of Theorem 4.1, until Section 6.4 



As in the conditions of Theorem |4.2[ suppose that 

M(x) > max{6 — x, x — a} 



for all x G [a, 6] and that m a = nib = 0. Then we may use Proposition 5.11 in place 
of Proposition |5.9| to the integrals in S3 and obtain 



(57) 
(58) 
(59) 



S, = - 



Htti 9(a)e(f(a) - ra) ^ y> g(b)e(f(b) - rb) 

^— ' O.irit f'(n\ — <i"l R— »rm ^— ' 



i?->oo z — ' 2ni(f'(a) — a 



^ 2m(f'(b) - r) 



+ lim V O f ^ 



C^(fl)/"(fl) 



+ 



U{a) 



+ 



'(a) -r| 3 M(a)(/ / (o) - r) 2 

r V n / U(b)f"(b) U{b) 
R ^r^( b ) Vl/'(fc) -H 3 M(6)(f (6) -r) 2 



We used a and 6 in place of a r and b r due to m a and m b being zero. Likewise, we 
have no a r , (3 r terms since by our assumption on M(x), a r = a and /3 r = b for all r. 
Just as in the proof of Theorem 4.1, the terms in line (57) equal 



-V(b) + V(a) + 0(A 2 (o) + A 2 (6)). 

we did the 
U(a)f"(a) 



The terms on lines (58) and (59) we bound as we did the terms in line (55), obtaining 

U(a) 



< 



M{a)\\f'{a)\\^ + ||/'(o)r 



17(6) 



+ 



t/(6)/"(6) 



*3 



M(b)\\f'(b)\\* 2 \\f'( 
< A 2 (a) + A 2 (6). 

Since we have no K r terms, we can let A 4 just equal 

U(x) ( rTTTi — \T.ri \\ ( l + \M'(x) 



f»(x)M(x) 



dx, 



and this completes the proof of Theorem 4.2 in this case. 
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8. Proof of Theorem 14.31 



The proof of this theorem mostly entails following the proof of Theorem |4.1| and 
simply being more careful with how certain error terms arise. 

we want to show that we can replace 0(T>% >1 (a)) + 0(Ai(a)) with 



In Section 



6.1 



some constant c\. But these terms simply serve as a bound for the sum 



00 pa 

Y I g(x)e(f(x) - rx) dx, 

„ J Or 



which, for fixed M(a), is a constant as b tends to 00. (If b is near a, we may need 
to take M(a) to be smaller, thus it's not always a constant.) Likewise, the term 
0(Ai(a) + A 2 (a)) that appear in Section 6.2 arise from bounding the term £(/3 r ) 



from Proposition |5.4| applied with the c and b of the statement of the proposition 
equal to a and (3 r , respectively. Provided b is large enough, (3 r will be constant, so 
the term £(/3 r ) may be likewise replaced by a constant c 2 . 

The error terms from Section |6.3| arise from bounding the sum 



(60) 



(frK) 



a<x r <b 



where £ r (a r ) (respectively, £ r (f3 r )) is the error from Proposition 5.4 applied with a 
(respectively, b) and c in the statement of the proposition equal to a r (respectively, 
(3 r ) and x r . 

We would like to replace the b in line (60) with 00. Consider first the £ r {ct r ) terms. 
Each such term is fixed (provided b is sufficiently large enough that b — a > M(a)). 
Thus we write 



E * 

a<x r <b 



a, 



\a<x r <b' b<x r <b'/ 

Suppose now b' is chosen to be on a sequence tending to infinity such that A3 (6') 
tends to 0. Then we have that 

lim > £ r (o: r ) 
b'^00 ^ 

a<x r <b' 

exists and is a constant c 3 by the assumed convergence of the integral in A 3 (a), and 
we have that 



b<x r <b' 



U{x) 



b f"{x){x-af 
V U(x) 



1 + 



1 



+ 



1 



dx 



+ 



l b f"(x)M(xf 
+ A 3 (6) + A 3 (fe') 
« A' 3 (b) + A' 3 (6') + A 5 



f"(x)M(x) f"(x)(x-a 

l + \M'(x 



1 + VFi^Mix] 



1 + 



f"(x)M(x) 



dx 
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by Proposition |5.14| Thus we have 

£ r (a r ) = c 3 + 0(A' 3 (b)) + 0(A 5 ). 

a<x r <b 

We must apply an additional idea to perform the same technique to the terms 
£ r ((3 r ). In particular, if x r + M(x r ) > b r , then (3 r changes as b grows. We may replace 
all such terms by £ r (x r ) with an error of size 0(A' A (b)). Following the argument of 
the previous paragraph, we have 



£r(Pr) = C 3 + 0(A' 4 (b))+0(A 5 ). 



a<x r <b 



In Section 6.4, the first sum on line (|46h is constant for b sufficiently large, and the 



second sum is just V°{b) + 0(A 2 (b)). 

By following the details of Proposition 5J3 that we applied to S3, we see that what 
we are left with at this point is 

hr(x) 



(61) 



(62) 



E 

f'(a)<r<f'(b) 



+ 



+ £ 

r<f'(a) or r>f'(b) J a 



27U 

br h r { X ) 
2tx% 



e(f(x) — rx) dx 
e(f(x)) — rx) dx 



where in the latter sum we assume we are taking the limit as R tends to infinity of 
the sum with the additional restriction that |r| < R. Let us, for ease of notation, 
rewrite the integrals as 



[a r ,b r ]\I Xr 

where as in condition (M) the set I Xr equals [x — M(x r ),x + M(x r )). We wish to 
replace b by a sequence of 6"s tending to infinity, such that A 2 (b') tends to 0. For 
such b' we have that above sum equals 

lim f ^2 ~ e(f(x) — rx) dx 



\ r \<R J \ <l r,K\Vxr 



lim > 

\r\<R 



[br,K]Vx r 



hr(x) 

2m 



e(f(x) — rx) dx. 



We now finish applying Proposition 5.9 by taking integration by parts of all the 



integrals and then bounding the variation. This combined with the ideas from earlier 



in the proof give us that the terms on (61) and (62) are bounded by 



c 4 + 0(A 2 (6)) + 0(A 5 ). 
This completes the proof of Theorem 4.3| 
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9. Proof of Theorem 14.41 



In the estimates of Theorem 4.1, we only roughly estimated the size of —T>(b) 



T>(a), which roughly are the first-order endpoint contributions at b and a, when /" 
was large. In this case, we included error terms of the size 0(U(a)) + 0(U(b)), which 
are at least the size of the first and last term of our initial sum. These terms can be 
improved, but require additional work. 



But recall, at the beginning of the proof of Theorem 4.1 we bounded a sum of the 
typd3 



(63) 



00 r-a 

/ g(x)e(f{x) — rx) dx 



in line (|27j) , which along with the sum in line (51) contributed to the 0(U(a)) error 



in V(a). In bounding this sum, we exploited the fact that the first-order endpoint 
contributions at a should cancel significantly (for example, in our estimation of line 



(29)). We expect, however, that more should be true, that not only should the 



first-order endpoint contributions display cancellations, but the integrals themselves 
should also display cancellation near a. 



As in the conditions of Theorem 4.4, assume that M(a), f"(a) > 1, let C and L be 
real numbers satisfying 

f"{a)- 1/2 < C < M{a) and Vf"(a) < L < f"(a) ■ min{l, C}, 

and let e := (a) and e' := (/'(a)). We now define a value a! r similarly to how we 
defined a r on line ( 26 ) . 



a + C if 1 < f'(a) -e'-r<L 
a-C if - 1 > f'(a) -e'-r>-L 
a otherwise 



If a' r 7^ a then a' r is between a and a r , except, possibly, for two values of r if L is very 
close to f"(a). 



We now break the original sum (63) into several pieces: 



E 

l<|/'(o)-e'-r[<L 



E 

l<\f'{a)-e'-r\<L 



E 

\f'(a)-e'-r\>L 



-0(A 1 (a)) 



where the last term comes from the term r = f'(a) — e' as in (28). In the sequel 



we shall refer to these three groups of integrals as Type I, Type II, and Type III 
integrals, respectively. 



11 



Since a r = a for almost all r, the sum is in fact finite. 
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We will postpone fully estimating the Type I integrals for now, and instead modify 
them into a more symmetric form. We first replace g(x) by 

g(a) + O I max g'(z) • \x — a\ J 

\ze[a' r ,a] J 

and fix) — rx by 

{f {a) - ra) + (f'(a) - r)(x - a) + \f\a){x - a) 2 



( max zeKia] /( 3 ^) 
+ O I L -^j J |x — a| 



Doing so gives 



/ 9(x)e(f(x) - rx) dx 

= g(a) £ e(f(x) - rx) dx + O £ \x - a\ ax) 

r 

= g{a) J e(/(x) — rx) dx + O 

J a' 



U(a)C 



M(a) 

9(a) £ e ((f(a) - ra) + (/'(a) - r)(x - a) + \f\a){x - a) 2 ) x 



x ( i+ Kif k - a|3,M * 

+ o ,u{a)c2 



M(a) 

= g(a) £ e ((/(a) - ra) + (/'(a) - r)(x - a) + ^/»(x - a) 2 ) dx 
j U{a)f»{a)C* \ 

v M( fl ) 

Shifting the bounds of integration and sending f'(a)—r to r + e', the Type I integrals 
can be rewritten as 

dx 



/a 
i g(x)e(f(x) - rx) 

g (a)e(f(a) + (e'-f'(a))a) ^ f 

1 r J a' —a 



" (.("(«) 2 , \ , 

e — - — x + rx + ra + e x dx 



+ 



l<\r\<L ' 

U(a)f"(a)C 4 r 
M(a) 
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The Type II and Type III integrals we evaluate using Proposition 5.11 This gives 



us 



(64) 



E 

,l<\f'(a)-e'-r\<L 

= E 

\f'(a)-e'-r\<L 

+ E 



+ E 

|/'(o)-e'-r|>Z,- Up 

g( a ' r ) 



2m(f>(a> r )-r] 
9(a) 



(65) 



(66) 



(67) 



2m(f'(a) -r) 

L<\f'(a)-e'-r\ U V 1 ' 

\f'(a)-r\<f"(a) 

o 



g(x)e(f(x) — rx) dx 

e (f( a 'r) ~ ra 'r) 

e(f(a) - ra) 



O 



+ 



E 

|/'(a)-e'-r|<L 

E 

L<\f'(a)-e'~r\ 
\[/'(a)-r|</"(o) 

E 

L<\f'(a)-e'-r\ 



U{a) 



U(a)f"(a) 

M(a)(f'(a' r )-ry \f'(a' r )-r\ s 

\ 



U(a) 



U(a) 



M(a)f"(a) M(aff"{af 



J 



U{a) 



U(a)f"(a) 

M(a){f'(a) -r) 2 ' |f(a)-r| 3 



+ 



Similar to Proposition 5.10 we have that 

f'(a±c) = f'(a)±f(a)-C(l + 



rj-C 
2 ■ M(a) 

with implicit constant 1, so that \ f'(a' r ) — r\ ^> f"(a)C for 1 < \ f'(a) — e' — r| < L. 
Thus, the sum in line (l64|) is bounded by 



E 

\f'(a)-e'~r\<L 



U{a) 



O 



f"(a)Cj \f"(a)C 



U{a)L 



and the sum in line (|65|) is bounded by 

U(a) 



° E 



,l<r 



U(a)f"(a) 

M(a)(f"(a)C + r) 2 (f"(a)C + r) 3 



O 

o 



U{a) 



+ 



U(a) 



M(a)f"(a)C f"(a)C 2 
U(a) 



()() 
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since c < M(a). The terms on line (66) are bounded by 



U(a) 



O 



U{a) 
M(a) 



M(a) f"(a)M(a) 3 , 

Finally, we apply Euler-Maclaurin summation to the sum on line (67) to bound it by 

'U(a)f"(a)' 



U(a) U(a)f"(a) 
M(a)L L 2 



O 



L 2 



since L is at most f"(a). 

Using all of these estimates for the integrals on line (63) (with the negative sign 
now, as they appeared in the proof of Theorem 4.1) and adding the terms from line 



(51), we obtain the following. 

oo 

- E 



E 



g(x)e(f(x) — rx) dx — hm 

\r\<R 
\f'(a)-r\>f"(a) 



g(a)e(f(a) - ra) 
2m{f'{a) - ra) 



-g(a)e(f(a) + (e'-f(a))a)x 



E 



(68) 



(69) 



i<M<£ gn(r) 

E 



x 2 + rx + ra + e'x ) dx 



> ' ? 9i ?l ^ 

r^oo Z— 2m(f'(a) - r) 

\r\<R ; 
L<\f'(a)-e'-r\ 



e(f(a) - ra) 



+ 



U{a)f"{a)C*L U{a)L U(a)f"(a) U{a) U{a) 



M(a) 



+ 



f"(a)c 



+ 



L 2 



f"(a)C 2 + M(a) 



Now we estimate the Type I integrals and the sum on line (|68j) in three separate 
cases. 

9.1. Case 1: a is an integer. In this case, we begin by adding the r and — r terms 
from the Type I integrals together, obtaining 



g(a)e(f(a) + (e'-f'(a))a) £ f 



7» 2 / \ , 

e | x + rx + ra + e x dx 



l<\r\<L ' 
"0 / til 



(70) 



—g(a)e(f(a)) 2_. \ e \ — 7, — %2 + rx \ ■ 2i sin(27re / x) dx. 

Kr<L V 2 / 



We estimate each integral in ( 70 ) by an application of integration by parts. 



° / /"(o) o 
( | x + rx ] • 2z sin(27rer) ax 



c 
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(71) 

(72) 

(73) 

The term 



sin(27re'x) 
n(f"(a)x + r) 



f"{o) 2 
x + rx 



c 



+ 



c 



2e' cos(2vre , x) f f"{a) 



o f"(a)x + r 



-x 2 + rx I dx 



C r/i 



f"(a) sin(2 7 re / x) ( f"(a) 
„ n(f"(a)x + ry € 

sin(27re';r) 



-a; + rx I fix 



ir(f"(a)x + r) 



on line (71) is when x = and bounded by 0(e' / f"(a)) if x = C. 

We will then apply the first and second derivative tests (Lemmas 5.6 and 5.7) to 
the integral on line (72). The total variation may be bounded by 

dx = O 



2e' cos(27re'x; 
f"(a)x + r 



e sin(27re'a;) 



O 



o 

+ 

r 



dx 



f"(a)x + r 
c e7"(a)cos(27re / x) 



o 



+ 



(f"(a)x + r) 2 
sin(27re'x)| dx 

f"{a) 
o (f"(a)x + ry 



dx 



dx 



O 



e'\(l + \e'\C) 



O 



Since | cos(27rex)| < 1 and f"(a)x+r is at least r for x G [0, C], the maximum modulus 
on the interval is at most \e'\/r. Therefore the integral on line (72) is bounded by 



O 



|e'|(l + |e'|C) 



mm 



We will also apply the first and second derivative tests to the integral on line (73). 
Using sin(27re'x) = 2ire'x + 0(e 3 x 3 ), we have that 

" c /» sin(27r e 'x) ff"(a) 



o n(f"(a)x + ry 

c 2e'f"(a)x 



o (f"(a)x + rf 



-x + rx \ dx 



/"(«) 



x + rx I dx + O 



C tn 



f"(a)e' 3 x 3 



dx 



(>2 
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Here we used that the function ax /(ax + 6) 2 is at x = 0, increases monotonically 
to a maximum at x = b/a, and then decreases mononotonically to as x tends to 
infinity. 

Thus the Type I integrals (70) are bounded by 
9(a) £ O 



Kr<L 



< \ V- ^ |e'|(l + |e'|C) f 1 1 
+ p(o) £ °('- 

Kr<L ^ 



7"(a)|e'| 3 C 4 



U(a)\e'\L 



O 



C/(a)|e'|(l + |e'|C)log(l + /"(«)) 



f"{a) 

+ 0{U(a)f"(a)\e'\ 3 C 4 ), 

which gives the bound on T>q(o) in this case. 
Also, we bound the sum in line (|68|) by 



O 



U{a) 



O 



U(a)f"(a) 
L 2 



using Proposition 5.16 



9.2. Case 2: f'(a) is an integer, a is close to an integer. As in case 1, we sum 
r and — r terms and use the first derivative test, obtaining 

g(a)e(f(a) + (e'-f(a))a)x 

'/"(a) 



x E 

l<\r\<L 



C-sgn(r) 



-x 2 + rx + ra + e'x ) dx 



/ /-// 



f"(a) 



f 

= 9(a)e(f(a) - f'(a)a) £ 2isin(27rer) / e 

l<r<L ^ C 

= (u{a) Hr-M 

V l<r<i / 

= 0(U(a) \e\L). 

By the same argument, we may complete the Fourier series in (68). 



x + rx I dx 



lim > 



^<l/'(a)-r| 



e(/(a) - ra) 
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g(a)e(f(a) - f'(a)a) e(er) 



27ii r-s>oo ^— ' r 

|r|<R 



■r<L 

= V(o)^(o)e(/(a) - f'(a)a) + O (l7(a)|e|L) . 
Together these give T>o(a) in this case. 

9.3. Case 3: f'{a) is an integer, a is far from an integer. In particular, by saying 
a is far from an integer, we want |e| > C . This implicitly requires that C < 1/2. 

In this case, we sum over all positive values of r seperately from all negative values 
of r. Ignoring the constant multiplier 

g(a)e(f(a) + (e - f'(a))a) 

for the moment, the sum over positive r values give 

"0 



Kr<L 



V- f (f"{°) 2 I \ J 

> / e x + rx + ra + e x ax 

f° (f'ip) 2 \ , , 

/ el — - — x I y, e{r{x + a)) ax 



l<r<L 





c 



(75) -/ c K^ + ( L + X ) ^nR^)) dX - 



Here we will apply the second derivative test to the integrals in lines (|74|)-( 75). 
The variation plus maximum modulus in both integrals is bounded by 



O 



1 



.(M-c). 

so that both integrals (as well as the corresponding integrals for the sum over negative 
r) are bounded by 

1 



O 



(|e|-C7)v7W 



In this case, again, we bound the terms in line (68) by 0(g(a)/L\e\). But since 
\J f"(a) <C L, this term is dominated by 



(N -c)/Fw 
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Thus we have T>q(o) in this final case. 

9.4. Combined case analysis. When e' = and M(a) < /"(a) 7 , the following 
choices of L and C optimize the error terms from line ( |69[ ): 

L = f"{afl l5 M{afl lb and C = /»- 2 / 5 M(a) 1/5 , 
when |e| < f , (a)~ 3/5 M(a)- 1 / 5 or /"(a)- 2 / 5 ]^) 1 ^ < | e | ; 

L = f"(a) 1/3 \e\- 1/3 and C = f"(a)- l \e\-\ 
when f / (a)" 3 / 5 M(a)- 1 / 5 < |e| < /' '{a)- 1 ' 2 - and, 

L = f"(a) 2/3 \e\ 1/3 and C = e/2, 

when f"{a)- 1 ' 2 < \e\ < f" (a)- 2 ^ M(afl\ 

Inserting these values into the error terms for cases 2 and 3 gives the result at the 
end of Theorem 14.41 

These methods work similarly at h in place of a. 

10. Proof of Corollary 11.91 

Recall that we are attempting to bound the error in the following van der Corput 
transform. 

(76) Yl* e ( (f ) = Y,* ^ ' e (~ 4r3 + X /8) + A 

n<N ^ / (l/12) 1 /2< r <(AT/i2)i/2 

We can apply the main theorems of this paper to either the sum on the left or 
the right of (76). However, we obtain best results if we use both at different times, 
dependent on how big || (N/12) lj/2 \\ is. 

10.1. Case 1: || (iV/12) 1 / 2 1| > (12iV) _1//4 or \\{N /Uf^W = 0. In this case, we apply 
Theorem 4.3 to the left-hand side of (76). Thus 

a = l, b — N, f( x )=(£j 3/ \ g{x) = l, 

G(x) = 0, H(x) = f^(x) = -^=x- 3 / 2 , W (x) = £- 2 , r (x) = 2^ 

and r' (x) is non-zero. For C 2 , C 2 ~, C 4 , D , D x , and D 2 equal to 2 and 5 equal to 1/2, 
we have that condition (M) is satisfied for M(x) — e ■ x for some small fixed e > 
(independent of N) and U(x) = 1. Hence H 2 — G > on J = [a(l — e), 6(1 + e)], so 
Jo is the full interval J, and both J± and J nu u are empty. 



We first need to check that the conditions of Theorem 14.31 hold. 

Using our assumption that ||(A^/12) 1 / 2 || > (12JV)- 1 / 4 , we can estimate A 2 (6) by 

A 2 (6) < A^- 2 (l + N 3 ^)(l + N- 1 ' 2 ) + N- 1 ' 2 
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^ + ^73 IIW12)'/1=0 

' 1 + 1 \\(Nll2) l / 2 \\> (\2NY l / A 

iV|| (iV/ 12) V2 1| 2 ^ ATi/2 1| (TV/12) V2 1|3 llv / > II ) 

< iV- 1/2 ||(iV/12) 1/2 ||*- 3 . 

Thus A 2 (6) tends to as N = 12k 2 and k tends to infinity along the positive integers. 
We can bound A3 (b) by 

< at- 3 / 2 , 

which clearly tends to as iV tends to infinity. 

The set K b equals [JV(1 + e) _1 , AT]. Let N be the largest integer in [1,N - 1/2] 
for which (AT/12) 1 / 2 is also an integer. Then the first integral and sum of A' 4 (6) are 
bounded by 



1 + x 1/2 + — dx 



[N(l+e)-\N] (N -X) 3 V JV-X 



AT JV(l + e) 

+ — — =- + 



-1 



(N-N) 3 (N — N(l + e) _1 ) 



3 



iV 1/2 N 



< — — =— + 



(JV - AT) 2 (AT - N) 3 
However, by the mean value theorem, we have, for some 77 G [N, N], 

N\ 1/2 fN\ 1/2 \ A f V \ 1/2 



N - N = [[ T2) "UJ J- 24 V12 

> \\{N/12) 1/2 \\* -N 1/2 , 

and A" x N, so that the first integral and sum of A 4 (6) is bounded by 
1 1 1 

ATI / O I I / TIT /-I /O I I ...O ATI /Oil / AT / -I \ 1 /OII...O 



A^l/2 1| (AT/12) 1/2 ||*2 ATl/2||(iV/12)l/2||*3 ~ N l/2 || ( iV / 12 ) 1/2 || *3 " 

Likewise the final integral and sum of A' 4 (6) are bounded by 



« f *- 5/2 (l + * 3/4 )(l + i) dx 

JNil+e)- 1 V X I J 



'N(l+e) 

+ AT" 2 (1 + AT 3 / 4 ) 
« AT" 3 / 4 . 
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Thus A' A (b) is bounded by 

< A^- 1/2 ||(iV/12) 1/2 ||*- 3 

and converges to on the sequence N = 12 • k 2 , k = 1, 2, . . . . 

Finally, we need to show that the integrals and sums in A3 (a), A4, and K, converge. 
The integrand of the integral in A3 (a) is 



x 



1 x 1 / 2 

1 + 



(x — l) 3 \ x 1 / 2 x — 1 
which converges like x~ 2 . The integrand of the integral in A4 is 

1 



x 



- 5 / 2 (l + x 3 / 4 ) 1 



1/2 



X 



which converges like x _7//4 . The integrand of the integral in K,(Jq, Wq^tq) is x~ 5 / 2 . 
The sum in the K term is bounded by |Wo(l — e)| + |Wo(-^(l + since Jo is a 
single interval and r' never changes sign. Since Wq(N(1 + e)) converges like N~ 2 , 
this shows that the conditions of Theorem |4.3| are satisfied. It also shows that A 5 is 
bounded by 0{N~^ 2 ). 

To apply the theorem we calculate T>(b) explicitly by 

'0(N- l l 2 ) || (TV/12) 1 / 2 1| = 

e ((iV/3) 3/2 ) x 



V(b) 



||(iV/12) 1 /2|| > (12iV)- 1/4 ' 

x '- MWi2m + *W<W 12 ^' 



If || (TV/12) 1 / 2 1| = 0, then A a (6) = N~ 2 ; and otherwise, we get that Ai(fo) = 0, since 
(127Y)- 1 / 4 is greater than f"(N) = (N/3)~ 1/2 /l2 for N > 1, and so m N equals 0. 
Combining the error terms together we see that 

A = 0(iV- 1/2 ||(7Y/12) 1/2 ||^ 3 ), 

which completes the proof of the corollary in this case. 



10.2. Case 2: 0^ || (TV/12) 1 / 2 1| < (12TV)- 1 / 4 . Now we apply Theorem|44|to e(-l/8) 
times the conjugate of the right-hand side of (76). We use the conjugate to guarantee 
that f"(x) is positive. 
We have 



IX. 



a = y/1/12, b = y/N/12, f(x) = 4x 3 , g(x) = VWa 
G(x) = -41472x, H(x) = UOVGx, H{xf - G(x) = 127872x. 
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Since G is never on [a, b] and H 2 — G is always positive on [a, b] we have that Jo is 
empty and J± = [a, b] . So, to complete our useful definitions, we have 



,/Q7 4- 7 

r ± (x) = 12(11 ± 2v / 37)x 2 , W ± (x) = Z. ' x~ 9/2 

and again r' ± (x) never equals on [a, b]. 

Again, for C2, C2-, C4, /Jo; Di, and /J2 equal to 2 and 5 equal to 1/2, we have that 
condition (M) is satisfied for M(x) = e ■ x for some small fixed e > independent of 
N and U(x) = g(x) = V2Ax. 

Since a is fixed, we have 

P(o)+Ai(a)+A 2 (o) = 0(l). 
By Theorem 4.4| we have that T>(b) equals 

- m 9 mm - r m + o ( M(6)a/ ^ (6)1/1B + §J) 
+o(f/(6)r(&) 1/3 ii&n 2/3 ) 

= -2^(v/iv7l2)(3Ar) 1/4 e(-(iV/3) 3/2 ) + 0(iV 3/2 °) 
+ 0(iV 5 / 12 ||(iV/12) 1 / 2 || 2 / 3 ). 
By our assumption that N is an integer, = and so 

Ai(6) < r 7/a < 1. 

For A 2 (6) note that in this case f'(b) is always an integer, and hence ||/'(6)|| equals 
0. Therefore, we have 

a 2 (6) « r 9 / 2 (i + 6 3 / 2 )(i + b) + r 3 / 2 + r 1 / 2 + r 1 / 2 
« iv- 1 / 4 

For A 3 (a) + A 3 (b) we may assume a = a + 1/2 and 6 = 6 — 1/2 (at worst, this 
assumption only makes the A 3 terms larger). Thus, 

,6 j 

' „-2 



A 3 (o) = [ — 7TT7 —r^ 1 + x 2 + — t— ■ 

KJ \J a+1/2 x^( x - 1/2)3 y x(x-l/2)> 

+ o(i) + o(r 9 / 2 ) 

= 0(1) 

and 

+ o(r 3 '' 2 ) + o(r 3 ) 



()8 
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3 (b-x) 3 {b-xY 

= 0(1). 

To estimate A4, note that by our above definitions we have 

\W-{x)\\r'_{x)\ + \W'_(x)\ + \W+(x)\\r'(x)\ + \W'(x)\ = 0(aT 7/2 ), 



and so, 



A 4 = o (^j\- 1/2 {l + x 3 / 2 )(l + x- 2 ) dx^j + x- 1 ' 2 dx^j 



+ |W_(a)| + |W_(6)| + \W+(a)\ + \W+{b)\ 

= 0(1) + 0(1) + 0(1) + o(b- g/2 ) + 0(1) + o(r 9 / 2 ) 

= 0(1). 

Thus we have shown that, in this case 

A = 2^( v ^v7l2)(3iV) 1 /4 e M ^ 3/2 + 0(iV 3 / 20 ) 

+ 0(iV 5 / 12 ||(iV/12) 1 / 2 || 2 / 3 ), 

by undoing our conjugation, multiplying by e(l/8), and moving this term to the 
correct side of the equation. 

11. Improving the results when ((/'(z 1 )!! is small 

In all the main results and the various examples given in this paper, the results 
tend to be at their best when ||/'(/x)|| (for \i equal to a or b) is large compared with 
\J f"(fi). However, when II/XaOH is on the order of y f"(fi), the estimates tend to 
be dominated by an inexplicit term, and even when ||/'(/x)|| is mu ch sm aller than 



a/ f"(/j), we still do not have great asymptotics, such as in Corollary |l .9 



Could the error term A in this range be made more explicit? If one's notion of 
explicit is a closed form that can be estimated moderately fast on a computer, then 
yes; if one's notion of explicit requires good asymptotic data, then the answer so far 
appears to be no. 

Suppose that f"(x) is generally smaller than 1 on [a, b]. The largest contributions 
to A come from bounding an integral of the form 

«n 

g(x)e(f(x) — rx) dx 



where \x r — /i| v//"(//). We will assume for a moment that x r is in the interval 
[a, 6]. 
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In the proof of Theorem 4.1, we bounded this integral with two techniques: first, 



possibly shifting fi to fi r = fi ± M (//) if \x r — fi\ was smaller than /"(/•*), and second, 



applying Propositions 5.4 and |5.5[ our variants of Redouaby and Sargos' results. 



However, another result of Redouaby and Sargos provides an alternate bound for 



such an integral. Under the conditions of Lemma 5A and the additional assumption 
that c—a <C N, M to simplify the error terms, we have the following result (combining 
Lemma 3 and Lemma 9 in [SO]) for a < c: 

g(x)e(f(x)) dx = -||=e(/( C ))J(^( fl ~c))+0 (i a) ~ 

vf ( c ) 

2 



„ (UM ( My /m 1/9 , xN 
+ °(^( 1+ atJ ( T y +( c -°)) + 




where a = /"(c)/2, 



/(c + f)-/(c) 



1/2 



= sgn(t) 

and denotes the modified Fresnel integral functio 

7{u) = / e(x 2 /2) dx. 
Jo 

An analogous result holds for the integral from c to b. 

Let us assume once again that M = N. Just as in Propositions 5.4 and 5.5. 
we see an error term of the form 0{UM/T 3 / 2 ). This new estimate on stationary 
phase integrals is useful — in the sense that the remaining error terms are of order 
of magnitude smaller than the explicit term — provided M ^> 1, T is of order much 
larger than 1 , and c — a is of order much smaller than 

M 1 ' 2 



mm 



f „ {c y/4 



This in turn is useful to the van der Corput transform if the good range for c — a 
includes c — a = /"(c)~ 1//2 , and this will be true provided f"(c) is of an order of 
magnitude much greater than M~ x . 

To apply these new estimates, we need to make a few further small changes in the 



proof of Theorem 4.1 First, for this r, we do not change /i to /x r , removing the term 
0(U(ji)/y/f"(ji)) from Ai(/i). Second, since we only apply the above estimates to 
an integral on one side of x r , we obtain terms of size [/(//) //"(/i)M(/i) = 0(A 2 (//)) 



from the terms now left uncancelled when we apply Propositions 5.4 and |5.5| to the 



other side. Finally, since we do not have any integral from /i to fx r in S3, all the 



15 More standard forms of the Fresnel integral functions would include a; 2 /4 or x 2 /2n where we 
have used x 2 /2 and would further break the function up into its real and imaginary parts. 
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associated terms no longer need to be bounded. If fi = a, then this means we lack 
the term 

g(a r )e(f(a r ) - ra r ) 
27ri(f'(a r ) — a) 

as well as 0(h r (a r )) and 0{h T (a r )). This thus removes the term with from 
T>([i), and reduces the last pair of terms in A 2 (/i) to 

When x r is not in the interval [a, b], say x r > b, then we need to be slightly more 
careful. We will need to assume that condition (M) holds on the larger interval [a, x r ], 
and then consider the difference of integrals 



-M(av) 



g(x)e(f(x)) dx. 



We would apply Propositions ^4 and 5J3 to the first integral and the above estimates 
to the second. 

While this may seem like a great improvement over earlier theorems — and in a 
sense, it is — it also is very hard to apply. The region ||/'(aOH x a/ /"(/•*) requires 
us to understand the function J-{u) in the region u x 1. Unfortunately, while there 
are very good asymptotics for J-{u) when u tends to or infinity, there are no good 
asymptotics when uxl. 

12. Further questions 



In Section 6J3, we do not make use of the fact that the intervals taken as arguments 
in the K r not only fail to include the point x r but also are excluded from an interval 
around x r . Is there a simple way to take this into account and does doing so improve 
the results in a substantial way? 

Can the results of this paper be explained in terms of the geometry of curves, in 
the same way as [1]? 

Dekking and Mendes- France conjecture in [5] that r(nlogn), the curve related to 
the exponential sum 

^2 e(nlogn), 

l<n<x 

is resolvable. Can the results of this paper be used to prove this? 
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